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In this review article we describe spin-dependent transport in materials with spin-orbit 
interaction of Rashba type. We mainly focus on semiconductor heterostructures, how¬ 
ever we consider topological insulators, graphene and hybrid structures involving super¬ 
conductors as well. We start from the Rashba Hamiltonian in a two dimensional electron 
gas and then describe transport properties of two- and quasi-one-dimensional systems. 
The problem of spin current generation and interference effects in mesoscopic devices 
is described in detail. We address also the role of Rashba interaction on localisation 
effects in lattices with nontrivial topology, as well as on the Ahronov-Casher effect in 
ring structures. A brief section, in the end, describes also some related topics including 
the spin-Hall effect, the transition from weak localisation to weak anti localisation and 
the physics of Majorana Fermions in hybrid heterostructures involving Rashba materials 
in the presence of superconductivity. 
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I. INTRODUCTION 

All the modern electronic devices are based on a clas¬ 
sical paradigm of negative carriers — electrons — and 
positive charge carriers — holes. This leaves completely 
unused an additional degree-of-freedom carried by both 
carriers: the spin. Spintronics explores phenomena that 
interlink the spin and the charge degree-of-freedom. It 
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is the field where traditional solid-state physics and ma¬ 
terial research have created their strongest bond, with 
each taking alternate leading roles. From our perspec¬ 
tive, spintronics can be mainly divided into two distinct 
subfields: metal and semiconductor spintronics (Zutic 
et ai, 2004; Sinova and Zutic, 2012). 

Metallic spintronic devices originate from the discov¬ 
ery of giant magnetoresistance (GMR) in 1988 (Baibich 
et ai, 1988; Binasch et ai, 1989) and the subsequent de¬ 
velopment of the spin valve (Dieny et ai, 1991). The 
GMR effect can be understood by assuming that any 
spin current is carried by two “types" of carriers, spin-up 
and spin-down. The two-channel picture of spin trans¬ 
port proposed by Mott explains the behaviour of magne¬ 
toresistive devices (Camley and Barnas, 1989), including 
GMR and tunnelling magnetoresistance (TMR) (Mood- 
era et al. , 1995), as well as spin injection into met¬ 
als (Johnson and Silsbee, 1988). 

In the GMR effect, two ferromagnetic layers sandwich 
a non-ferromagnetic metal layer of nanometer thickness. 
When the magnetisation of the two ferromagnetic layers 
is parallel, the valve is in a low resistance state (open). 
When the two are antiparallel, the valve is in a high re¬ 
sistance state (closed). The TMR effect rather than the 
GMR one, takes place when a metallic spacer is substi¬ 
tuted by an insulating barrier. In TMR at room tem¬ 
perature, the spin-vale effect increases by a factor of ten 
with respect to GMR. This is the basic principle of hard 
disk read heads and magnetic random access memories. 
These are some of the most successful technologies of the 
past decades, with scaling trends outdoing even comple¬ 
mentary metal-oxide-semiconductor, the technology for 
constructing integrated circuits. Albert Fert and Peter 
Griinberg have been awarded the Physics Nobel Price in 
2007 for their studies on the GMR effect. 

Semiconductor spintronics, on the other hand, has 
the potential to be integrated seamlessly with nowadays’ 
semiconductor electronics. It is progressing along a sim¬ 
ilar path as metallic spintronics and has achieved a re¬ 
markable success in the past decade. The spin-field ef¬ 
fect transistor (spin-FET) was proposed in 1990 (Datta 
and Das, 1990) as a development to the GMR set-up 
where the two ferromagnets are left collinear and the “on" 
and “off" status is achieved by rotating the injected spin 
when travelling through the two contacts. This geometry 
does not require a magnetic field for switching the rela¬ 
tive magnetisation of the two ferromagnets, thus reduc¬ 
ing spurious magnetic fields into the electronic circuits. 
The spin-FET is an electronic analog of an electro-optical 
modulator for photons. A ferromagnet injects spin po¬ 
larised carriers into a semiconductor channel. Here, due 
to the modulation of the spin-precession length, it is pos¬ 
sible to vary the spin orientation. This modulated spin 
signal is then detected by a second ferromagnet. The 
spin precession length is varied via top- and back-gate 
voltages that tune the strength of the spin-orbit inter¬ 


action (SOI) of Rashba (R) type (Rashba, 1960; Yu A. 
Bychkov and Rashba, 1984). The RSOI originates from 
the lack of inversion symmetry along the growth direc¬ 
tion of the semiconductor heterostructure that is hosting 
a two-dimensional electron gas (2DEG) (de Andrada e 
Silva et al., 1994, 1997). It can be tuned by changing 
the shape of the confining potential via the application 
of an external electric field (Nitta et al., 1997; Schapers 
et al., 1998; Grundler, 2000; Miller et ai, 2003; Meier 
et al., 2007). However, the conductivity mismatch prob¬ 
lem (Johnson and Silsbee, 1988; Schmidt et al., 2000) 
between a ferromagnet and a low-dimensional semicon¬ 
ductor has hindered, for a long time, an efficient injec¬ 
tion of spin polarised currents. Progress for solving this 
problem came with the research on ferromagnetic semi¬ 
conductors (Ohno, 1998). Nevertheless, only recently it 
has been possible to implement completely the function¬ 
ality of the spin-FET first in Ref. (Koo et al, 2009), and 
after, with minor changes, in Ref. (Chuang et al., 2014). 

Semiconductor spintronics is also relevant because it 
has permitted to observe other coherent phenomena, 
such as the Aharonov-Casher (AC) effect (Aharonov and 
Casher, 1984) in ring-type and other loop structures (Yau 
et al., 2002) so as proposed by many groups (Nitta 
et ai, 1999; Frustaglia and Richter, 2004; Bercioux et ai, 
2005a; Capozza et ai, 2005; Lucignano et ai, 2007). The 
Aharonov-Casher effect is the analog of the Aharonov- 
Bohm (AB) effect (Aharonov and Bohm, 1959) for parti¬ 
cles with an angular momentum in an external electrical 
field (Richter, 2012). In the case of carriers in a semi¬ 
conductor 2DEG, the RSOI permits the coupling of the 
electron spin with the external electric field. Other phe¬ 
nomena, e.g., the spin-Hall effect (D’yakonov and PereF, 
1971; Hirsch, 1999; Murakami et ai, 2003; Sinova et al, 
2004), have been predicted and observed in semiconduc¬ 
tors (Kato et al, 2004; Sih et ai, 2005a; Wunderlich 
et ai, 2005a, b; Sih et ai, 2006a) before their observa¬ 
tion in metals (Valenzuela and Tinkham, 2006). A key 
role in the success of semiconductor spintronics is the 
tunability of the RSOI via external gating. 

The material in this review article is presented in the 
following form: in Sec. II we present a summary of the 
materials with sizable RSOI effects, these range from the 
standard semiconductor heterostructures, carbon-based 
materials, topological insulators (TIs) to Weyl semimet¬ 
als (WSs). In Sec. Ill, in order to introduce some basic 
transport property of carriers in systems with RSOI, we 
revise the spectral properties of 2DEGs and graphene 
with RSOI, in addition we spend some effort to under¬ 
line the basic difference in the presence of a further lat¬ 
eral confinement, thus moving from the two-dimensional 
(2D) to the quasi-one-dimensional-case (quasi-ID). Af¬ 
ter, we introduce the concept of spin-double refraction 
that is a phenomenon taking place in a hybrid struc¬ 
ture composed of region with RSOI sandwiched between 
two regions without RSOI. In Sec. IV we give a defini- 
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tion of pure spin currents and propose a derivation of 
a formula for a spin current in the standard Landauer- 
Biittiker approach, furthermore we introduce two pro¬ 
posals for generating pure spin currents — based on the 
ratchet and the pumping effect. Next, we give an in- 
depth analysis of the spin-FET by Datta and Das and an 
explanation of its working principle, finally we present a 
recent experiment showing a possible realisation of such 
a transistor. In Sec. V we give a complementary look 
at the effects induced by RSOI, specifically we will anal¬ 
yse its effects as non-Abelian gauge held. Particularly, 
we look at non-Abelian effects in quantum networks and 
quantum rings. In Sec. VI we give a very brief and gen¬ 
eral overview over other aspects of the physics associated 
to the RSOI, specifically we look at the spin-Hall effect 
(SHE) in Sec. VI.A, than, in Sec. VI.B, we will see how 
RSOI can modify the weak localisation (WL) into weak 
anti-localisation (WAL) in a disordered system. Finally 
in Sec. VI.C we will quickly review the vibrant held of the 
quest of Majorana quasi particles in quantum wires with 
RSOI, magnetic held and superconductivity. In Sec. VII 
we present our outlook for this growing held of condensed 
matter physics. 


II. MATERIALS WITH RASHBA SPIN-ORBIT 
INTERACTION 

On a very general basis, SOI routes down to the rela¬ 
tivistic correction to hydrogen-like atoms (Bransden and 
Joachain, 2003; Sakurai and Napolitano, 2014). It reads 

Uso = - 2H s-( E (r) xp) (1) 

[ZjTTLC ) 


where p is the momentum operator and s = (s x , s y , s z ) 
is a vector of Pauli matrices describing the spin S = 
operator. By considering a Coulomb electrostatic poten¬ 
tial V(r) = — ^r- we have that the strength of this rel- 
ativistic correction goes as thus proportional to the 
atomic number Z (Bransden and Joachain, 2003; Saku¬ 
rai and Napolitano, 2014). From the relativistic correc¬ 
tion expressed in Hamiltonian (1) we learn that materials 
characterised by a higher atomic number will present in 
general a stronger SOI, this will be the case for TIs and 
Weyl- and Dirac-semimetals. 


A. Semiconductors heterostructures 

In this review we shall mostly discuss the physics of 
low-dimensional semiconductor heterostructures. The 
simplest example is a 2DEG. This is commonly realised 
in III-V semiconductor heterostructures by modulating 


the doping density along the growth direction — usu¬ 
ally (001) (Bastard, 1988; Seeger, 2004). In solids, SOI 
can have different nature depending on the crystal sym¬ 
metries (Samokhin, 2009), however in this review arti¬ 
cle we shall mostly deal with the two main indepen¬ 
dent sources of SOI in III-V semiconductors, namely the 
RSOI, due to the lack of structural inversion symmetry 
(SIA) (de Andrada e Silva et al., 1994, 1997) and a Dres- 
selhaus (D) term due to a lack of bulk inversion symmetry 
(BIA) (Dresselhaus, 1955). The three-dimensional (3D) 
SOI correction to the free electronic Hamiltonian due to 
BIA reads: 

^DSOI Px (Py ~ Pz) s x + Py ( P z ~ Px) s y 

+ Pz (pl ~ P 2 y ) S z . (2) 

In order to achieve the effective Hamiltonian acting on 
the electrons confined in the 2DEG we have to inte¬ 
grate along the growth direction with the constraint that 
(p z ) = 0 , while (p 2 z ) jz 0 is a sample dependent constant. 
Therefore we have 

%DSOI = 0 (PySy - PxSx) + 7 (p x p 2 yS x - p y pls y ) (3) 

where 7 is a material dependent constant, and /? depends 
on (pi). The first term contains the so called linear DSOI 
whereas the second one describes the cubic one (Winkler, 
2003). Usually the latter can be disregarded as (pi) -C 
Px ? Py i n a 2DEG. 

Due to charging effects, the quantum well confining 
electrons in two dimensions is never perfectly flat, there¬ 
fore at the interface between the two differently doped 
semiconductors, a potential gradient arises W = E. 
This effective electric field couples with the electrons mo¬ 
tion as 

■Hrsoi oc (E x p) ■ s . (4) 

Considering a quantum well along the (001) growth di¬ 
rection E = E z z, the RSOI is rephrased as: 

a ex 

^RSOI = X P) ■ Z = -(p y S x -p x Sy) , (5) 

where a depends on the material and on the confin¬ 
ing potential. By simple inspection one can notice that 
Rashba Hamiltonian (5) and linear Dresselhaus SOI (3) 
are equivalent under the unitary transformation rotating 
the spin Pauli matrices s x —> s y , s y —► s x and s z —> —s z 
therefore the spectral properties of electrons subject to 
RSOI or linear DSOI are exactly the same (Lucignano 
et ai, 2008). It is important to notice that through a 
simple mathematical manipulation we can recast RSOI 
in the form s ■ -B e ff (p) where B c g oc (z x p) is an effective 
magnetic field rotating the carrier spin in the plane of the 
2DEG. 

Despite this formal equivalence, the role to RSOI has 
been investigated with much more attention in the last 
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FIG. 1 (a) Schubnikov-de Haas oscillations as function of the 
gate voltages (Nitta et al, 1997); (b) ted conduction band 
diagram and electron distribution. Panels (a) and (b) cour¬ 
tesy of Ref. [Nitta et al, 1997]; (c) Experimental magneto¬ 
conductance, Act = a(b) — ct(0) (circles), offset for clarity, 
along with three-parameter fits (solid line) for several gate 
voltages. Inset: Experimental magneto-conductance data for 
the most negative gate voltage, showing pure weak localiza¬ 
tion. (d) Density and mobility as function of , extracted 
from longitudinal and Hall voltage measurements; (e) Exper¬ 
imental conductivity, showing strong dependence on V g . Pan¬ 
els (c)-(e) courtesy of Ref. [Miller et al., 2003]. 


decades for two main reasons. First, in conventional III- 
V 2DEGs a > /3. Second, more interestingly, from few 
decades a can be easily tuned by gating the heterostruc¬ 
ture (de Andrada e Silva et al., 1994, 1997) to some hun¬ 
dreds of meV nm. The spin precession length can be 
expressed as 

7T TtH 2 

ko 1 * ? 

kso am* 

in III-V quantum wells it ranges from hundreds of nm 
to few /i?7i, thus comparable with typical dimensions 
achievable in the fabrication of modern mesoscopic de¬ 
vices. The RSOI was first predicted in bulk semicon¬ 
ductors by Rashba (Rashba, 1960). Later Bychkov and 
Rashba (Yu A. Bychkov and Rashba, 1984) proposed to 
use the Hamiltonian (5) in order to explain some mag¬ 
netic properties of semiconductor heterostructures (Stein 
et al., 1983; Stormer et al., 1983). In figure 1(b) we 
show the shape of the quantum well and of the elec¬ 
tron density profile showing the position of the 2DEG. 
The lack of inversion symmetry in semiconducting 2DEG 
heterostructures is at the origin of the RSOI. Therefore a 
change in the asymmetry can also lead to a change in the 
strength of the RSOI. For example, by applying a gate 
voltage is possible to modify the profile of the potential 
confining the 2DEG with a secondary effect of changing 
the strength of the RSOI. Thus, in a first approxima¬ 
tion the coupling constant a depends on the gate voltage 
a —> a(V). 

A change in the strength of the a, thus, can produce 
a change in the beating pattern of the Schubnikov-de 


Haas oscillations (de Andrada e Silva et al., 1994, 1997) 
[c.f. Figure 1(a)]. This effect was first observed by the 
group of Nitta (Nitta et al., 1997) and after by other 
groups (Schapers et al., 1998; Grundler, 2000). An¬ 
other method for probing the variation of the coupling 
constant is to observe the transition from WL to weak 
WAL (Miller et al., 2003; Akkermans and Montambaux, 
2007) 1 , see for instance figure l(c)-(e). 


B. Quantum wires 

By applying a further confinement to a 2DEG it is pos¬ 
sible to confine the electrons in a quasi-ID structure — 
a quantum wire. A detailed investigation of the RSOI in 
InGaAs/InP quantum wires of different width has been 
put forward by Schapers et al. [Schapers et al., 2004, 
2006]. Quantum wires can also be made out of InAs, 
GaAs, InSb or other materials (Nilsson et al., 2009; Nadj- 
Perge et al., 2012; Estevez Hernandez et al., 2010). Re¬ 
cently, quantum wires with sizable RSOI are attracting 
a lot of interest as they are the fundamental building 
block for many proposals aiming at the observation of 
Majorana quasi particles in condensed matter physics. 
In these setups the interplay between RSOI, supercon¬ 
ductivity and magnetic fields (Oreg et al., 2010; Lutchyn 
et al., 2010) is pivotal. 2 

C. Carbon-based materials 



FIG. 2 (Colours online) Spin-resolved photoemission spectra 
along TK in the vicinity of the Fermi energy, (a) Overview 
spectra and (b) 7r states at arbitrary offset. Courtesy of 
Ref.[ Varykhalov et al, 2008]. 


1 More details about this phenomenon will be given in Sec. VI.B. 

2 Some more analysis of Majorana physics is reported in Sec. VI.C. 
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The RSOI has been predicted also for single 
layer graphene (SLG), using simple symmetry argu¬ 
ments (Kane and Mele, 2005a, b). A more rigorous 
calculation has been proposed by Huertas-Hernando et 
al. (Huertas-Hernando et al, 2006), where a modulation 
via atomic Stark effect and curvature effects were con¬ 
sidered. However, due to the small atomic number of 
the carbon atoms [see equation (1)] the strength of the 
coupling a is very small in the SLG structure. There 
are proposal for overcoming this obstacle and increasing 
the strength of the SOI, by coating the graphene surface 
with ad-atoms thus inducing a spin-dependent hopping 
mediated by the surface impurities. This is the case, for 
example, of SLG deposited on a Ni substrate and in¬ 
tercalated with Au atoms. Spin-resolved angle resolved 
photoemission spectroscopy measurements have shown a 
sizable RSOI of the order of some meV (Varykhalov et al., 
2008; Dedkov et al, 2008; Li et al., 2011; Marchenko 
et al., 2012). Similar proposal have been advanced also 
for hydrogen impurities (Castro Neto and Guinea, 2009) . 
The modification to the energy spectrum of SLG due 
to RSOI are shown in figure 2. Large RSOI has been 
also proposed by using rotating magnetic fields (Klino- 
vaja and Loss, 2013). 

The RSOI in SLG is described by a Hamiltonian very 
similar to (5), the main difference is that it depends only 
on the pseudo-spin and not on the momentum operators. 
To the lowest order in the momentum expansion: 

^!o G = \ ( s x > (6) 

where s are the Pauli matrices describing the electron 
spin and <x its pseudo-spin. Here the pseudo-spin is the 
internal degree of freedom describing the presence of two 
inequivalent carbon atoms in the honeycomb lattice of 
SLG (Neto et al, 2009). 

The RSOI has been predicted also in bi-layer graphene, 
due to the interaction with a substrate of and an external 
electric field (Gelderen and Smith, 2010; Guinea, 2010; 
Konschuh et al, 2012; Mireles and Schliemann, 2012), 
however so far it is has not been measured. 

In order to conclude this overview of carbon-based ma¬ 
terials with RSOI, relevant effects are predicted also for 
carbon nanotubes (CNTs) (Martino et al, 2002; Martino 
and Egger, 2005; Izumida et al, 2009). In Refs. [Kuem- 
meth et al, 2008 and Jespersen et al, 2011] the RSOI 
interaction strength is measured in CNT quantum dots. 

D. Topological insulators 

Topological insulators are band insulators hosting spin 
polarised edges states within the bulk gap (Hasan and 
Kane, 2010; Qi and Zhang, 2011; Ando, 2013). 

The original idea about this new state of matter was 
put forward by Kane and Mele in two seminal research 


papers published in 2005 (Kane and Mele, 2005a, b). 
They showed that a specific type of SOI in SLG can lead 
to the opening of a gap in the spectrum and that edge 
states exists within this gap when considering a ribbon 
geometry. In figure 3 we can observe the gap opened 
by the intrinsic SOI (Huertas-Hernando et al, 2006) and 
the edge states within the gap. In two-dimensions, this 
phenomenon is also known as quantum spin-Hall effect 
(QSHE). It can be thought as two copies (one per spin) 
of the integer quantum Hall effect (IQHE). The main dif¬ 
ference is that, contrary to the magnetic field, SOI does 
not break time-reserval symmetry. Thus, here, there are 
two counter propagating helical edge states opposed to 
the single chiral edge state of the quantum Hall fluid. As 
in the case of the IQHE, the edge states are characteris¬ 
tics only of the finite size systems. It is also possible to 
introduce a specific topological numbers named Z 2 that 
defines the robustness of these edge states (Kane and 
Mele, 2005a). 

Up to now, such phenomenon has not been observed 
in SLG, however shortly after the seminal works of 
Kane and Mele, a proposal for observing the QSHE 
in HgTe/CdHgTe quantum wells was put forward by 
Bernevig et al. [Bernevig et al, 2006 and Bernevig and 
Zhang, 2006] , whose prediction was readily confirmed by 
Konig et al [Konig et al, 2006]. Since than, there is 
a huge quest for discovering new 2D TIs. One of the 
most interesting is represented by InAs/GaSb quantum 
well (Liu et al, 2008; Knez et al, 2011) which is very at¬ 
tractive as it is achievable using the standard and, very 
well known, III-V semiconductor technology. 

Contrary to the case of IQHE, the paradigm of the 
QSHE can be extended to 3D. This work has been mainly 
carried out by Fu and Kane [Fu et al, 2007 and Fu and 
Kane, 2007]. Since their seminal works, a race for discov¬ 
ering new 3D TIs is open. Recent overviews on these ma¬ 
terials can be found in one of the several review articles 
that has been published on the topic (Hasan and Kane, 
2010; Qi and Zhang, 2011; Ando, 2013). Among them 
the paper by Ando (Ando, 2013) contains an interesting 
table summarising all known 2D and 3D TI discovered 
up to 2013. 

Most of the research work is devoted to analyse the 
spectral properties of such material—mainly via angle- 
resolved photoemission spectroscopy. Only few experi¬ 
ments show the transport properties of this new class of 
materials. Mainly because most of the 3D TI are not 
proper insulators due to a non negligible current flowing 
through impurity states (Bardarson and Moore, 2013). 
A very interesting set of transport experiments has been 
performed in films of BUSea (Qu et al, 2011) prepared 
with a special two-dimensional geometry (c.f. figure 4). 
In this work the authors show a sizable tunability of 
RSOI in 3D TI. The effects of this tunability are clearly 
observed in the interference pattern of the Aharonov- 
Bohm (AB), the Altshuler-Aronov-Spivak (AAS) and the 
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FIG. 3 (Colours online) Edge states of the quantum spin-Hall 
effect in zig-zag graphene nano ribbons, the bulk states are 
given by the solid-orange lines, the spin-polarised states are 
the dashed-blue lines. 


Aharonov-Casher (AC) effects (Akkermans and Montam- 
baux, 2007; Altland and D., 2010; Nazarov and Blanter, 
2013). This is very interesting because it would allow to 
realise novel functionalities such as possible non-Abelian 
operations on spins. 


E. Weyl semimetals 

Weyl semimetals also known as “topological semimet¬ 
als” are zero gap semiconductors with a Fermi level very 
close to the centre of the gap that is at zero energy as 
in the case of SLG and TIs. The effective Hamiltonian 
describing the low energy states close to the Fermi level 
in WSs is represented by a 2 x 2 Hamiltonian that is 
linear in the momentum — as for TI and SLG — but 
it is an actual 3D model, thus containing all the three 
Pauli matrices. The simplest model Hamiltonian in kp 
approximation reads: 

T~L± = ±vfp ■ s , (7) 

where uf is the Fermi velocity, and p and s are the vectors 
of the momentum and of the Pauli matrices, respectively. 
By inspection of equation 7 we evince that a generic per¬ 
turbation cannot gap the energy spectrum. Its energy 
spectrum reads E±(p ) = ±up|p|. This is usually correct 
for an even number of nodal points where the expansion 
(7) is allowed. At each node we can associate a chiral¬ 
ity, this measures the relative handedness of the three 
momenta and the Pauli matrices associated in the Weyl 
equation — the chirality ±1 for the Hamiltonians T~L±, 
can also be thought as a source of Berry curvature car¬ 
ried by the WS nodes. This is a general property of Weyl 
fermions realised in band structures: their net chirality 



FIG. 4 (Colours online) (a): SEM image of the square-ring 
device realised in BLSes; (b): Illustration of the AB effect in¬ 
terference (red trajectories) and AAS interference (solid and 
dotted green loops) for charges in a square ring. The exis¬ 
tence of RSOI creates an effective magnetic field Bso pointing 
towards/ or from the center of the ring for counterclockwise 
(CCW)/clockwise (CW) propagation modes, which influences 
spin precession and generates an AC phase in addition to the 
AB and AAS phases, (c) Ai( 8 ) and A 2(8 ) as a function of 
8. Here 8 = 2am*L /ft 2 is spin precession angle over a dis¬ 
tance L and A\ ( 8 ) = (cos 4 8 + 4 cos 8 sin 2 8 + cos 28) /4 and 
A 2 ( 8 ) = (sin 2 8 + cos20)/4 are the modulation of the wave 
function as a function of 8 (c.f. Ref. (Qu et al, 2011) for more 
details). The former varies roughly at twice frequency of the 
latter, (d): A R — V g curves at fixed magnetic fields marked 
by the lines and arrows of corresponding colours. The spin 
precession angle is modulated by 47r by varying V g from the 
interval 2.16 to 2.77 V marked by the two stars. The dashed 
lines help illustrating the opposite phases between the green 
and blue curves in the AAS region. Courtesy of Ref. [Qu 
et al, 2011 ]. 


must in fact cancel 3 . Note that it is usually assumed that 
bands are individually nondegenerate. This requires that 
either the time-reversal symmetry or the inversion sym¬ 
metry (parity) is broken. In order to realise the minimal 
case of just a pair of opposite chirality Weyl nodes, time- 
reversal symmetry must be broken (Wan et al., 2011). 

The idea of WS was first theoretically proposed by Mu¬ 
rakami in 2007 [Murakami, 2007] and later it was further 
elaborated by Wan et al. in 2011 [Wan et al, 2011]. 
Based on aft-initio band calculations, it was proposed 
that pyrochlore iridates, such as Y 2 H 2 O 7 , in the anti¬ 
ferromagnetic phase, may realise such a WS [Wan et al . , 
2011 and Yang and Kim, 2010]. A recent review on the 
spectral properties of this type of system is reported by 
Vafek and Vishwanath [Vafek and Vishwanath, 2014] and 
on the transport properties by Hosur and Qi [Hosur and 
Qi, 2013], 


3 This can be understood in terms of the fermion doubling theo¬ 

rem (Nielsen and Ninomiya, 1981). 
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III. SPIN TRANSPORT IN RSOI MATERIAL 

Before describing the role of RSOI in the spin- 
dependent transport, it is important to understand how 
it modifies the spectral properties of a free electron in a 
2DEG. 

A. General properties of the RSOI in semiconductors 

We consider a 2DEG in the (x, y)-plane in the presence 
of the RSOI (5) and with effective electron mass to *. The 
full Hamiltonian of the system is 

•p 2 a , 

w » = w + (8) 


The eigenvalues read 

fc2 

E±{k)= 2nS ±ak= 2 ^ (fc±fcso) Aso ’ (9) 


A; 2 is the modulus of the electron 


where k = \k\ = ^jk 2 

momentum, kso = is the RSOI coupling constant 

with the dimension of a momentum and Ago = x • 

Usually the last term of (9) is neglected because it is sec¬ 
ond order in a. Moreover, even if included, it would lead 
to a rigid shift of the bands, thus renormalising the chem¬ 
ical potential p i. The eigenvectors of the Hamiltonian (8) 
relative to the spectrum (9) are plane waves: 



FIG. 5 (Colours online) Properties of the Rashba energy spec¬ 
trum. (a) Portion of the energy spectrum of the Hamiltonian 
(8). (b) The Fermi contours relative to the Hamiltonian (8), 
the spin states are shown as well, (c) Section of the energy 
spectrum for a free electron, (d) Section of the energy spec¬ 
trum for an electron in presence of a magnetic field, e.g. Zee- 
man splitting, (e) Section of the energy spectrum for an elec¬ 
tron in presence of RSOI. 


velocities are given by 

V±(k) = 1 = -*r ± ak = ± kso)k. (11) 

h ok to* to* 

If we consider the quantum-mechanical velocity operator 


T ± (r) 


e ifc r / 1 \ 

V2 ' v ±ie- i V 


( 10 ) 



( 12 ) 


where 9 = arctan(fcj / /fc a; ) is the polar angle between the 
momentum vector and the k x direction. It is important 
to note that the spin states (10) are always perpendicu¬ 
lar to the motion direction. In fact, if an electron moves 
along x direction the spinor part of the eigenvectors be¬ 
come (1; i) and (1; —i) that is the spin up and spin down 
are along the y direction. By contrast, if the electron 
moves along the y direction, the eigenvectors become 
(1; 1) and (1; — 1) that is the spin up and spin down state 
in the x direction [see figure 5(b)], 

In figure 5(c)-(e) we report the fc y -section of the energy 
spectrum vs the momentum for a 2DEG in different phys¬ 
ical situations. In 5(c) we considered a free particle in a 
2DEG. In this case the spectrum is twofold degenerate 
in spin. In the presence of a magnetic field B [see panel 
5(d)], the spin degeneracy is lifted out by the Zeeman 
effect and the gap separating spin up and spin down is 
equal to 2g*y,BB, where g* is the effective gyromagnetic 
ratio and is the Bohr’s magneton. When the RSOI 
is present 5(e), the spin degeneracy is lifted out except 
for k y = 0. In this situation the degeneracy is removed 
without the opening of gaps. The semiclassical particle 


and the expression (10) for the eigenstates, it is straight¬ 
forward to show that its matrix elements are given by 

($±(fc)| r |'F ± (fc')) = S(k - k') v±(k). (13) 

That is the semiclassical velocities v±(k) are, as usual, 
the diagonal elements of the velocity operator. In the 
presence of the RSOI, the velocity operator is not sim¬ 
ply the momentum divided by the effective mass as for 
free electrons and the Fermi surface splits into two sur¬ 
faces shown in figure 5(b). One can easily observe that 
counter propagating electrons have opposite spins. Spin 
and momentum are locked to each other. This will have 
importance consequences in the following. Parametrising 
the wave vectors as k = k (cos p, sin tp) the two Fermi 
surfaces are described by the following equation: 


fc±(£ F ) = =F 


a to * 
H 2 



(14) 


Here the double sign corresponds to the two dispersion 
branches (9). 




















B. Rashba SOI in quantum wires: subbands hybridisation 


120 


0.5 


By further confining a 2DEG along one direction (we 
choose y in the following) we can realise a quantum 
wire. In this case we cannot solve analytically the sys¬ 
tem Hamiltonian, as we have done for the simple case 
of a 2DEG. Different theoretical models have been used 
to describe confinement effects. A very convenient rep¬ 
resentation consists in assuming a transversal confining 
potential in the y directions, and let the electrons propa¬ 
gate along the the x direction. Such choice is particularly 
convenient when expressing the matrix elements of the 
SOI (Moroz and Barnes, 1999; Governale and Ziilicke, 
2002; Perroni et ai, 2007). The Hamiltonian reads: 


H = H ll +H ± +H mix , (15) 


with the following terms: 


Hu = 

n ± = 

^mix 


Px 

2m* 

P 2 y 

2m* 

a 


hkso 

m* 

+ V(y) 


&xPy ■ 


(16a) 

(16b) 

(16c) 


Where V (y) is a infinite well potential of width W. Ne¬ 
glecting "Hmix, the terms in equations (16a) and (16b) do 
commute, therefore the eigenvalues and eigenvectors can 
be easily expressed as: 


-®||+_L (^b kx) 


h 2 / 2 7r 2 ?V 

2 m* V X + ~W 2 


tpna y kAx,V) = <t>n(y)lpk x ( x )Wv) 


h 2 ksok x 

m* 


Wy) 

(17) 

(18) 


where 4> n {y ) are the ro-th eigenfunctions of the potential 
V ( y ) that are a either a sine or a cosine functions and 
ipkx (x) simple plane waves shifted by the RSOI, here ( a y ) 
can get the values ±1. The term ?{ m i x induces mixing 
between these states and gives rise to a deformation of 
the electronic bands and to anti-crossings in the energy 
spectrum. 

Anti-crossings occur between sub-bands corresponding 
to transverse eigenstates with different n [figure 6 (a)]. 
While simple one dimensional (ID) models predict a rigid 
spin-momentum locking with in-plane spin perpendicu¬ 
lar to the momentum along the wire direction, multi¬ 
band models predict that only electrons with momen¬ 
tum far away from anti-crossings have the spin essen¬ 
tially perpendicular to the momentum, i.e. (S y ) ~ 
and k ~ k x x. This is well shown in figure 6 (b), where 
we can see the expectation value of (S v ) as a function 
of the longitudinal momentum k x — it changes form — 1 
to 1 for k x changing from negative to positive values. 
The same picture applies to the other sub-bands. Here 
we shall focus on the role played by the number of sub¬ 
bands used for evaluating the system properties. The 





-10 -50 5 10 
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FIG. 6 (Colours online) (a) Band structure of quantum wires 
with square well potential along the confining direction. The 
coloured lines refer to the first three sub-bands in the case of 
a wire with three modes whereas the dotted-black lines refer 
to the same wire but only considering two modes, (b) Spin 
polarisation (S y ) as a function of the transversal momentum 
for the first sub-band in the case of a two band model (red- 
solid and red-dashed lines), (c) Spin polarisation (S y ) as a 
function of the transversal momentum for the second sub¬ 
band in the case of a two band model (red-solid and red- 
dashed lines). Both in panel (b) and (c) the spin polarisation 
is evaluated using a truncated Hilbert space containing 50 
modes (blue-solid and blue-dashed lines). 


Hamiltonian (15) cannot be diagonalized exactly, but a 
partial analytical/numerical solution can be obtained by 
truncating the Hilbert space associated to the modes pro¬ 
duced by V(y). Diagonalizing the Hamiltonian (15) for 
N 1 modes produces a result that is substantially cor¬ 
rect for the lowest N — 1 modes but not for the mode 
N. This is observed, e.g., in the spectral properties and 
in the spin polarisation. The spectrum is shown in fig¬ 
ure 6(a), we have in colours-solid lines for N = 3 and 
in black-dashed lines for N = 2. We see that the first 
sub-bands in both cases (N = 3 and N = 2) are quite 
similar apart from a negative energy shift (Perroni et al . , 
2007). However, the situation is very different when we 
consider the second sub-bands, these greatly differ in the 
two cases — due to the presence of the third mode for 
the case N = 3. Considering the spin polarisation the 
effect is more visible. In figure 6(b) the spin polarisa¬ 
tions for the first sub-bands for N = 2 (red lines) and 
N = 50 (blue lines) are shown, the behaviours are simi¬ 
lar in both cases but around k x ~ 0. The error produced 
by Hilbert space truncation is shown in 6(c), here it is 
shown the spin polarisation for the second sub-bands (red 
lines), for N = 2. We observe a mirroring of the results 
reported in 6(b). However, the most correct result with 
N = 50 (blue lines) contradicts completely N = 2. 

The message here is that when solving a problem of 
quantum transport in a confined geometry a sufficient 
number of sub-bands should be taken into account in or¬ 
der to avoid systematic errors due to the truncation of 
the Hilbert space (Moroz and Barnes, 1999; Governale 
and Ziilicke, 2002; Perroni et ai, 2007). The same con- 
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elusions are obtained by changing the shape of the con¬ 
fining potential V ( y ) from square potential one to har¬ 
monic oscillator one (Moroz and Barnes, 1999; Governale 
and Ziilicke, 2002) or by introducing a further periodic 
modulation along the longitudinal direction x (Smirnov 
et ai, 2007). The sub-band hybridisation that we have 
shown before can give rise to dip into the conductance 
of quantum wires with RSOI as shown in several works 
both analytically (Moroz and Barnes, 1999; Perroni et al . , 
2007) and numerically (Sanchez and Serra, 2006) and can 
also be interpreted in terms of a Landau-Zener transition 
when analysing the conductance of a quantum point con¬ 
tact (Eto Mikio et ai, 2005). 

C. General properties of the Rashba SOI in graphene 

If we restrict ourself to the case of a single valley ap¬ 
proximation, the effects of RSOI on SLG via Hamiltonian 
(6) change with respect to what we have seen in the case 
of a 2DEG in Sec. III.A. The complete Hamiltonian for 
the SLG with RSOI reads: 

H k = Vpcr ■ p + - (s X <j) z 

OL 

= V F (<?xPx + VyPy) + ^ ( S ^ tT V ~ S x a y) ■ ( 19 ) 

The energy eigenstates are plane waves r ip ~ $(fc)e lfc ' T ’ 
with $ a four-component spinor and eigenvalues given 
by (v F = h= 1) 

£±,e(k) = + fcy + —, (20) 

where index e = ± specifies the particle/hole branches of 
the spectrum. The energy dispersion, as a function of k x , 
at fixed k y = 0, is illustrated in figure 7 with and without 
RSOI. As we can see here, the RSOI gives rise to a finite 



FIG. 7 (Colours online) Energy spectrum of graphene for 
k y = 0 as a function of k along the x axis. The dashed lines 
correspond to the linear dispersion for a = 0 , the solid and 
the dotted lines to the case of finite RSOI. The lines with the 
same colours correspond to the same spin state. 

curvature of the linear energy spectrum and lifts the spin 
degeneracy. Two opposite spin bands get gapped while 
the other two are still degenerate at k ~ 0. However, 
contrary to the case of 2DEGs, the energy spectrum is 
not shifted along the momentum axis as in equation (9). 
This is related to the fact that the RSOI Hamiltonian 
does not depend on the electron momentum but only on 
its pseudo-spin, to the lowest order in the momentum 
expansion. Nevertheless, this approximation no longer 
holds if we expand RSOI Hamiltonian to higher orders, as 
they explicitly depend on the electron momentum. Their 
main effect, for strong RSOI, is to produce the appearing 
of extra Dirac cones (tridiagonal wrapping Ref. [Zarea 
and Sandler, 2009; Rakyta et ai, 2010; and Lenz et ai, 
2013]). 


The eigenstates of the Hamiltonian (19), are expressed 
by the spinors d>± je (fe): 


*±, e (*0 


1 

2-^/cosh Q± 


x (e e0± / 2 , ee £ 0 ± / 2 ,±iee e6l± / 2 , iie 11 ^ eS ±/ 2 ), 


where D T denotes transposition and 

sinh 6 » ± =±^, ( 21 ) 

gi<fr _ k x -\-\k y ( 22 ) 


with k = yjk\ + ky. The spin operator components 
are expressed as Sj = (g> I 2 . Their expectation values 


on the eigenstates <I > ±, (r read 


,o \ esin< / ) 

1 x) 2 cosh 9± ’ 

(23a) 

e cos (f> 

{ v> ~ 2 cosh 0 ±’ 

(23b) 

(S x ) = 0 , 

(23c) 


which shows that the product ±e coincides with the sign 
of the expectation value of the spin projection along the 
in-plane direction perpendicular to the direction of prop¬ 
agation. For vanishing RSOI, the eigenstates e reduce 
to linear combinations of eigenstates of S z . Similarly, the 












expectation value of the pseudo-spin operator er is given 

by 
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_ \ _ ECOS <j> 

> x > cosh 6± 5 

(24a) 

' \ e sin 4> 

~ a yl ~ cosh 9± ‘ 

(24b) 


Since the RSOI in SLG does not depend on momentum, 
the velocity operator still coincides with the pseudo-spin 
operator: v = r = i [H,r] = cr. Thus, the velocity ex¬ 
pectation value in the state is given by equations 
(24a) and (24b). Alternatively, it can be obtained from 
the energy dispersion as 

ek 

v e = V k £ e = . (25) 

V fc2 + T 

If we considered also the intrinsic SOI, the effective Fermi 
velocity would acquire a more complex dependence from 
the RSOI (Lenz et ai, 2013; Bercioux and De Martino, 
2010 ). 

Sub-band hybridisation due to RSOI also occurs in 
graphene nano-ribbons. However, its manifestation de¬ 
pends on the type of boundary (Neto et al. , 2009): zig¬ 
zag (Zarea and Sandler, 2009) or armchair (Lenz et al . , 
2013). The most peculiar effects are observed in the case 
latter. Here in fact, the RSOI couples not only different 
modes but introduces also a finite coupling inside each 
sub-band. It is due to the presence of two inequivalent 
carbon atoms in the unit cell of a SLG nano ribbon and 
is completely different from the case of a quantum wire 
with RSOI of the previous section (Lenz et ai, 2013). 


D. The spin double refraction 

The physics of an interface between a 2DEG with and 
without SOI is very similar to that of the optical bire¬ 
fringence (Born and Wolf, 1999) and is a direct conse¬ 
quence of the two possible Fermi velocities for the two 
modes (11). In the following we consider a 2DEG in the 
(a;, j/)-plane characterised by an effective electron mass 
to* and an interface along the y direction separating a 
region without SOI (N-region) and a region with it (SO 
region). Elastic scattering at the interface allows for con¬ 
servation of the total energy and of the momentum par¬ 
allel to the interface — k y in our case. As we have seen in 
the previous section, the Fermi surface for an energy £p is 
constituted of two circles with radius (14). An incoming 
particle from the N region is characterised by the mo¬ 
mentum k N and the incidence angle </> = arctan 
Momentum conservation implies that 

sin 9+ = sin 9- = k N sin <f >, (26) 

where we have introduced k^P = fc^(fp) and 6± are the 
propagation angles in the SO region for the two different 




FIG. 8 (Colours online) Spin Refraction angles 9± as a func¬ 
tion of the incoming particle angle 4>. Panel (a) and (b) case 
parameters of two modes always open. Panel (c) and (d) case 
parameters where the critical angle of mode (-), is bigger 
than 7 t/ 2. In Panels (b) and (c) the red-dashed line is the 
mode (-) whereas the blue-solid line is the mode (+). 


modes. An incoming particle from the normal region 
propagates into the SO region along two different modes 
with indices ±. In Refs. [Perroni et al. , 2007; Ramaglia 
et ai, 2003; Khodas et al. , 2004; and Ramaglia et al . , 
2004] this phenomenon is named spin-double refraction. 
From equation (26) we can obtain an expression for the 
two refraction angles: 


/ fc N \ 

6± = arcsin ^^gQ-sin^J . (27) 

If the refraction angle is equal to |, than the correspond¬ 
ing mode will not propagate forward, thus it is closed. Of 
course, there are two critical angles <j>± corresponding to 
the closure of the corresponding modes in the SO region 
defined by the relation: 

/ /..SCI \ 

<t>t = arcsin ■ ( 2g ) 

According to equation (14): so that </>+ < <f>~. 

Furthermore, considering that is always smaller than 
fc N , the mode (+) will be always open. The previous re¬ 
sults are based only on kinematic considerations and are 
independent by the nature of the SOL The key ingredient 
is the presence of an interaction that splits the 2D Fermi 
surface in two circles. For example, the same physics 
would be possible at an interface between a normal re¬ 
gion and a region with a strong Zeeman splitting, assum¬ 
ing that the Fermi energy allows the propagation along 
the two modes in the region with the magnetic field. 

Some care is required when writing down the Hamilto¬ 
nian for this scattering problem. Indeed the two regions 
(N and SO) are obtained by considering an inhomoge¬ 
neous coupling constant for the RSOI a(x). As a con¬ 
sequence, a(x) does not commute with the momentum 
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operator [a(x),p x ] ^ 0. Thus the Hamiltonian describing 
the two regions has to read: 


■%n-SO = 7 ^— + 


a(x) 


2m' 


h 

1 da(x) 
’ y 2rri 


(SxPy SyPx ) 

+ Vi n tS(x). 


dx 


(29) 


where the second to last terms reestablish the hermitic- 
ity (Messiah, 2014). In the simplest approximation we 
can consider a(x) = aQ(x — rc; nt ), where 0 (x) is the 
Heaviside step function. The last term mimics the 
presence an SOI interface, according to the standard 
approach describing interfaces with delta potentials of 
strength V]nt (Blonder et aL, 1982). A more realistic de¬ 
scription has to account also for possible changes of the 
effective electron mass m* —► m*(x) = m^Q(xi n t — x) + 
m SO®( a; int + x ) i n the two regions, more details can be 
found in Ref. [Perroni et aL, 2007 and Ramaglia et aL, 

2003] . When considering the presence of a second inter¬ 
face, where the RSOI vanishes again, the refraction an¬ 
gles for the two modes will be the same and equal to the 
incoming one <p (Khodas et aL, 2004; Ramaglia et aL, 

2004) . However this double interface structure cannot 
produce a spin polarisation when the incoming particles 
are unpolarised. This is due to the fact that in order to 
produce a spin polarisation we need to violate the On- 
sager relation for a two terminal system (Adagideli et aL, 
2012; Gorini et aL, 2012). Violation is possible if and only 
if time reversal symmetry is broken, this implies, e.g., the 
application of a magnetic field. 

The same physics is possible also in the case of SLG, 
the main difference is that the RSOI in this case does not 
contain the momentum operator ( 6 ). As a consequence 
we do not need the extra term in the full Hamiltonian 
as in equation (29). A complete analysis of the spin- 
double refraction for the case of graphene is reported in 
Ref. (Bercioux and De Martino, 2010) where a transfer 
matrix method is presented, which is useful for study¬ 
ing the spin-dependent transport in hybrid structures in 
SLG, as a periodic RSOI potential (Lenz and Bercioux, 
2011; Costa et aL, 2013). 


IV. PURE SPIN CURRENT GENERATION IN RSOI 
MATERIAL 

A very crucial point in spintronics is the creation of 
pure spin currents. Over the last 20 years, many meth¬ 
ods for creating pure spin currents have been proposed. 
We will focus mainly on ratchet and pumping methods 
that have received also experimental verification. Here 
we define what is a spin current and discuss a method 
for evaluating it in terms of the Landauer-Biittiker for¬ 
malisms (Scheid et aL, 2007a). A pure spin current is 
defined as a particle flow carrying finite spin polarisation 
without an associated charge current. 


A simple example follows: Suppose, for instance, that a 
charge Q and spin polarisation S pq moves from a contact 
p to q during the time t £ [0, T] and in the next period 
t £ [T, 2 T] the same amount of charge Q, but a different 
spin polarisation S qp , moves from contact q to p. The 
net charge transported between p and q, in the time t £ 
[0,2T] is zero, but the net spin polarisation is AS = 
S pq S qp . 

In order to evaluate the spin-current we consider N non 
ferromagnetic contacts injecting spin-unpolarised current 
into the leads. We use, as customary, a local coordinate 
system for the lead under investigation, where x is the 
coordinate along the lead in the direction of charge prop¬ 
agation due to an applied bias in linear response and y 
is the transverse coordinate. In each lead, at a fixed en¬ 
ergy E, several conducting modes are open. The wave 
function for each mode reads 

y) = -TT^^^XniyMs) , (30) 

V kn(E ) 

here we have introduced the transverse eigenfunctions of 
the lead of width W: 

Xniy) = ^lw siu ( r w) VneN 

with the eigenenergy E n = (yp-) 2 and £(s) is the 

spin eigenfunction. The superscript ± of 4' refers to 
the motion direction along the lead axis with the wave- 
vector k n = \j2m*[E — E n )/h 2 . We use the scattering 
approach for deriving the current formula: the ampli¬ 
tudes of the states inside the leads are related via the 
scattering matrix S(E), determined by the Hamiltonian 
of the coherent conductor. The scattering state inside 
the q- th lead reads 

V>e( x : y) = Yl [ a ls( E )'&E,ns( X ’ 2 /) + b ls( E )^E,ns( X ^ V) ■■ 


with (s = ±). It consists of incoming states 4' + entering 
the coherent conductor from the contact q and outgoing 
states 4 ,_ leaving the coherent conductor into the contact 
q. The amplitudes of incoming a J ns and outgoing b l ns 
states are related each other via the scattering relation 

N 

K’A E ) = E E E Si Js>,ns( E ) a L( E )> ( 31 ) 

j= 1 n(Ej 5==tl 


where the scattering matrix S(E) has the following struc¬ 
ture for an N terminal system: 


S(E) 


r^\E) 

t l ' 2 {E) ■■ 

■ t^ N {E) 

t 2 ’\E) 

r 2 ’ 2 (E) ■■ 

■ t 2 ’ N {E) 

t N,1 (E) 

t N ’ 2 (E) ■■ 

■ r N ’ N {E ) 
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Here the sub-matrix r J,J (E) is a square matrix of dimen¬ 
sionality M J (E), corresponding to the number of open 
channels in lead j at energy E connected to a reservoir 
with chemical potential Hj — it already includes the 
spin degree of freedom. The matrix r 3 ' 3 (E) contains 
the scattering amplitudes of incoming channels of lead j 
being reflected back into outgoing channels of the same 
lead. The sub-matrix t l,3 (E) is a M l (E ) x M 3 (E) matrix 
that contains the scattering amplitudes for transmission 
between incoming channels from lead j and outgoing 
channels of lead i. 

In order to evaluate the spin-current we consider N non 
ferromagnetic contacts injecting spin-unpolarised current 
into the leads. We use, as customary, a local coordinate 
system for the lead under investigation, where x is the 
coordinate along the lead in the direction of charge prop¬ 
agation due to an applied bias in linear response and y 
is the transverse coordinate. In each lead, at a fixed en¬ 
ergy E, several conducting modes are open. The wave 
function for each mode reads 

y) = - 7 f=e ±ik ^ x Xn(yMs) , ( 32 ) 

V k n( E ) 

here we have introduced the transverse eigenfunctions of 
the lead of width W: 

Xn(y) = /J s m(^) Vn£ N 

with the eigenenergy E n = (yp-) 2 and £(s) is the 

spin eigenfunction. The superscript ± of U/ refers to the 
motion direction along the lead axis with the wave-vector 

k n = yJ^k(E — E n ). We use the scattering approach 
for deriving the current formula: the amplitudes of the 
states inside the leads are related via the scattering ma¬ 
trix S(E), determined by the Hamiltonian of the coher¬ 
ent conductor. The scattering state inside the g-th lead 
reads 

<Pe( X ’ y) = Yl [ a la( E )^E,ns( X ^ V) + b ns ( E )^E,ns ( X ’ V) ■ 
ns 

with (s = ±). It consists of incoming states 'I' + entering 
the coherent conductor from the contact q and outgoing 
states leaving the coherent conductor into the contact 
q. The amplitudes of incoming a J ns and outgoing b l ns 
states are related each other via the scattering relation 

N 

K'AE) = E E E Ws', ns ( E )aL(E), ( 33 ) 

j—1 n£j s==tl 

where the scattering matrix S(E) has the following struc¬ 
ture for an N terminal system: 

/r 1 ’ 1 ^) t 1 ’ 2 ^) ••• t 1 ’ N {E)\ 
t 2,1 (E) r 2 ’ 2 (E) ••• t 2 ’ N (E) 

\t N ’\E) t N ’ 2 (E) ■■■ r N ’ N {E)/ 


Here the sub-matrix r ( E ) is a square matrix of dimen¬ 
sionality M J (E), corresponding to the number of open 
channels in lead j at energy E connected to a reservoir 
with chemical potential /j ? — it already includes the spin 
degree of freedom. The matrix r 3 ’ 3 (E) contains the scat¬ 
tering amplitudes of incoming channels of lead j being 
reflected back into outgoing channels of the same lead. 
The sub-matrix t l ’ 3 (E) is a M l (E) x M 3 {E) matrix that 
contains the scattering amplitudes for transmission be¬ 
tween incoming channels from lead j and outgoing chan¬ 
nels of lead i. 

The wave function of the scattering state inside lead i, 
where only the incoming channel (ns) £ j is populated 

(ai's' = n<>s', s )’ reads for 3 = 

A E ,ns( X iy) =^E,ns( X ^) 

+ r n's',ns( E )^E,n's'( X ^y)> ( 34 ) 

(n / s / )Gi 


and, correspondingly, for ) 

VE,ns( x iy) = tl n's',ns( E )^E,n'A X ^y)- ( 35 ) 

(n's')Ei 


For a generic spin wave function $(a :,y) the spin current 
/g(x) passing a cross section (x = const) of a lead is given 
by: 


h( x ) = J dy^*(x,y)J s ^(x : y). (36) 


Here, we have used the most common definition of the 
spin current operator (Rashba, 2006). This is defined 
with respect to an arbitrary quantisation axis u, than 
the operator Js reads: 


Js 



h 

( 9 

_2m*\ 

\dx 



(37) 


In this operator we have the spin operator projected 
along the u quantisation direction times that standard 
quantum mechanical current. Here, as customary, the 
partial derivatives act on the expressions to their right 
and left (as indicated by the arrows), respectively. By 
acting this the former operator on the scattering state 
(34) we obtain for the spin current (36) inside lead i nrT 


na( x 6 *) = 


h 2 

2m* 


(n' 


n' s' ,ns 


(E) 


where (n, s £ j, j = i). For the other scattering state (35) 
we find the corresponding expression (n, s £ j,j ^ i) 


J S ;E,ns 


(x £ i) 


h 2 


* 


E 

(n's'jGi 


n' s' ,ns 


(E) 


S(E) = 


2m 


s 
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Each channel is populated according to the Fermi-Dirac 
distribution f(E\ fi n ) of the respective contact n, the 
statical average of the total spin current in lead i reads 


<,!( ' €i, ) = s?i dE 


N 


E E f( E '^o) I i,E,ns( X ^ i ) 

i =1 (ns) ej 




(38) 
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where the spin resolved transmission and reflection are 
defined as 
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Probability conservation implies that the scattering ma¬ 
trix has to be unitary [S(E)JS(E) = S(E)[S(E)J = I, 
than the following relation holds: 
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By using the symmetry relation between the spin trans¬ 
mission and reflection (c.f. Ref. [Scheid et ai, 2007a]), it 
is straightforward to show that 


R s i (E) + ^2A q (E) = 0 . 

q^i 


In view of Eq. (38) we eventually find for the spin current 
in lead i 

1 f°° 

h(x€i) =—I d EY'lfiE^iAfiE^qKTAiE). 

0 q^i 

(41) 

Equilibrium spin currents can locally exist in systems 
with SOI as shown for 2DEGs (Rashba, 2003) and in 
mesoscopic systems (Nikolic et ai, 2006), however equa¬ 
tion (41) clearly shows that at thermal equilibrium (fj,j = 
/ 1 i Vj £ {1,N}) the spin current vanishes inside leads 
without SOI and/or magnetic field. This absence of equi¬ 
librium spin currents in the leads has been shown for sys¬ 
tems with preserved time-reversal symmetry (Kiselev and 


Kim, 2005). An expression very similar to (41) has been 
used to investigate the SHE (Hankiewicz et al., 2004; 
Ren et al., 2006; Bardarson et ai, 2007; Briine et al . , 
2012). In the presence of SOI or magnetic fields it is 
nontrivial to write down a continuity equation for the 
spin current, the problem has been addressed in several 
articles (Nikolic et al., 2006; Sun and Xie, 2005; Shi et ai, 
2006; An et ai, 2012). 


A. Pure spin current generators 

In this section we describe two mechanisms that can 
be used for generating a pure spin current. 


1. Spin Ratchet 

A particle ratchet (Reimann, 2002; Hanggi and March- 
esoni, 2009) is a system with broken inversion (left/right) 
symmetry that generates a net currents upon external 
ac-driving in the absence of a net (time-averaged) bias 
potential. Ratchets have much in common with cur¬ 
rent rectifiers, though there are differences, in particular 
in the dissipative case (Reimann, 2002). The theoreti¬ 
cal concept of ratchets, originally introduced for classi¬ 
cal dynamics, was later extended to the quantum dis¬ 
sipative regime (Reimann et al., 1997). The main dif¬ 
ference between ratchets and rectifiers is that quantum 
ratchets exhibit current reversal upon changing, e.g ., the 
temperature or energy (Hanggi and Marchesoni, 2009; 
Reimann et al., 1997). Such quantum ratchets were ex¬ 
perimentally realised in semiconductor heterostructures 
in a chain of asymmetric ballistic electron cavities in the 
low-temperature regime, with a dynamics was close to 
coherent (Linke et ai, 1999). Also the charge current re¬ 
versal phenomenon has been demonstrated (Linke et al, 
2002 ). 

A spin ratchet is a quantum ratchet with lack of bro¬ 
ken spatial symmetry — in order to get a zero charge 
current — and a finite RSOI ensuring a breaking of the 
spin symmetry thus allowing for a finite spin current. 
Theory predicts spin ratchets to work both in the ballis¬ 
tic regime (Scheid et al, 2007b, 2010; Ang et ai, 2015) 
and in the dissipative one (Smirnov et ai, 2008a, b, 2009). 
A set-up for the ballistic ratchet has been also proposed 
by using periodically disposed magnetic stripes instead 
of electrostatic barriers and RSOI (Scheid et ai, 2007a, 
2006). 

A ballistic spin ratchet is mainly constituted by a 
quasi-one dimensional channel with a set of symmetric 
barriers and RSOI [see figure 9(a)]. The lack of bro¬ 
ken symmetry implies the absence of a net charge cur¬ 
rent. In order to understand the presence of a finite spin 
current we show here the same argument proposed in 
Refs. [Scheid et al., 2007b, 2010]. A key ingredient is 
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the sub-band hybridisation that we have introduced in 
Sec. III.A. 

We consider a wire with two open transverse modes 
(N = 1,2) and a smooth symmetric potential barrier 
U(x) in the two bias or rocking situations (V ^ 0), see fig¬ 
ure 9(b). Upon adiabatically traversing the barrier from 



ABC 



FIG. 9 (Colours online) Panel (a): Sketch of the system of 
symmetric periodic equally spaced barriers used for the bal¬ 
listic spin ratchet. Panel (b): Illustration of the spin polarisa¬ 
tion mechanism for transmission through a strip with a single 
adiabatic symmetric potential barrier U ( x ) (solid line) in the 
two rocking situations (dashed and dotted line). At points 
A, B and C the position-dependent energy dispersion relation 
E n (k x ) is sketched with respect to the Fermi energy Up (hor¬ 
izontal line) for two transverse modes and RSOI-induced spin 
splitting of each mode. Courtesy of Ref. [Scheid et al, 2007b] 


region A via B to C [see lower part of figure 9(b)], the 
electron energy spectrum E n [k x ) split by RSOI is shifted 
rigidly upward and downward by the adiabatic potential 
barrier U[x). For fixed Fermi energy Ep, the initial shift 
causes a depopulation of the upper levels (N = 2 ) and 
a spin-dependent repopulation of the lower one (N = 1) 
while moving from B to C. If Ep is traversing an anti¬ 
crossing between successive modes (see the region indi¬ 
cated by the dashed window in figure 9(b), there is a 
certain probability P for electrons to change their spin 
state. This causes an asymmetry between spin-up and 
spin-down states for the repopulated levels (Eto Mikio 
et al., 2005). The related transition probability can be 
computed in a Landau-Zener framework. In the case of 
a transverse parabolic confinement of frequency ujq, this 
probability reads (Eto Mikio et al., 2005; Scheid et al., 
2007b): 


P(±U 0 ) = l-exp 


— 7r/csoWoS 


-1 


d x [U(x,z) ± U 0 g(x,z)\ 


(42) 


Here E, denotes the difference in the polarisation of the 
two modes involved, and the function g(x,z) describes 
how the potential drops inside the system. In principle 
g{x, z) is obtained by solving the Boltzmann equation, 
however its exact form is irrelevant in order to understand 
the working principle of the spin ratchet. The spin trans¬ 
mission is proportional to P(V) and thus different for the 
two rocking situations V ^ 0. Hence, the ratchet spin 
current (/s(V)) is finite, also for the case of a symmetric 
barrier. A quantitative explanation of the spin ratchet 
effect for a system of multiple barriers [figure 9(a)] is be¬ 
yond this model. An experimental verification of the spin 



FIG. 10 (Colours online) The spin-field effect transistor pro¬ 
posed by Datta and Das [Datta and Das, 1990]. 


ratchet effect was proposed in Ref. [Costache and Valen¬ 
zuela, 2010]. However, in this experiment the breaking 
of the spin symmetry was achieved by combining a su¬ 
perconductor with a magnetic field. 


2. Spin Pumping 

Adiabatic charge pumping (Brouwer, 1998) consists of 
the transport of charge obtained — at zero bias voltage 
— through the periodic modulation of some parameters 
in the scattering region, e.g. some voltages. 

In 2003, Governale et al. [Governale et al., 2003] pro¬ 
posed an updated of the scheme of charge pumping ex¬ 
tended to the spin, following a suggestion by Mucciolo 
et al. [Mucciolo et al., 2002], The pumping scheme con¬ 
sists of an electrostatic barrier that is changed periodi¬ 
cally in time and a gate, as second pumping parameter, 
to tune the strength of RSOI (Governale et al., 2003). 
They showed that, neglecting sub-band hybridization, 
the charge current is zero whereas the spin current is 
finite. However, in a more realistic case in which many 
modes are opened, and consequently coupled by RSOI, 
the charge current is finite but two order of magnitude 
smaller than the spin one. A prototypal spin pump¬ 
ing system has been realised in a GaAs quantum dot in 
Ref. [Watson et al., 2003], in which an ingenious method 
for detecting spin currents is presented. Similar pump¬ 
ing mechanisms have been proposed for graphene with 
RSOI (Bercioux et al, 2012), and TIs (Citro et al, 2011; 
Inhofer and Bercioux, 2013; Hofer and Biittiker, 2013; 
Hofer et al, 2014). 


B. Datta and Das spin Field Effect Transistor 

The spin momentum locking due to the RSOI is one 
of the key ingredients used in 1990 for the proposal of 
a spin-FET by Datta and Das [Datta and Das, 1990]. 
It is believed that spin-FET has the advantages of low 
energy consumption and fast switching speed since it 
does not involve creating or eliminating the electrical 
conducting channels during the switching, required by 
traditional FETs. This spin-FET is the electronic analog 
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of an electro-optic modulator. The electro-optical effect 
makes the dielectric constant of a medium different along 
the two perpendicular directions ( e.g ., y,z). 

Let us assume to have photons with a polarisation at 
45° with respect to the y axis (in the y-z plane) (Born 
and Wolf, 1999). They can can be represented as a linear 
combination of 2 - and y-polarised photons: 



Finally, the analogue of the electro-optic material is 
realised by employing a ID quantum wire with RSOI. 
Within a single band model, RSOI causes the | f) and 
14-) electrons with the same energy to have different wave 
vectors and fcj. [see figure 5(e)]. Let us consider an 
electron travelling in the x direction with k z = 0 and 
k x ^ 0. The eigenenergies for the corresponding ID case 
are given equation 9 with k y = 0. 

Using equations (9)-(14), it is possible to recover a re¬ 
lation for the phase shift between the two spins: 


Because of the anisotropy of the dielectric constant, 
as the light passes through the electro-optic mate¬ 
rial of length L, the two components acquire different 
phase shifts k\L and kiL. Thus the light emerging 
from the electro-optic material can be represented as 
(e' klL , e lk2L ) . It is collected by an analyser polarised 
along the 45° direction. Therefore the output power is 
given by 


P 0 (X 
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= 4 cos 2 
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which is proportional to a. 

The above analysis was originally limited to a single¬ 
mode quantum wire. However, as we have learnt in the 


k= o t so * o 4 = o 
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The light output is modulated with a gate voltage that 
controls the differential phase shift Ad = {k\ — k 2 )L. 

The analog electronic device based on the RSOI is 
shown in figure 10. In the original proposal by Datta and 
Das, the polariser and the analyser were suggested to be 
implemented via ferromagnetic contacts (as Fe) (John¬ 
son and Silsbee, 1988). In such materials the density of 
states for electrons with a specific spin orientation — at 
the Fermi energy — greatly exceeds that for the opposite 
direction. A contact magnetised in the x direction prefer¬ 
entially injects and detects electrons spin polarised along 
positive x which is represented as a linear combination 
of the positive z-polarised | f) and negative ^-polarised 
electrons |j,): 



FIG. 11 Schematic of the tight-binding model for the system. 
In the shaded areas the spin-orbit interaction if finite tso ^ 0. 
Courtesy of Ref. [Mireles and Kirczenow, 2001]. 


Sec. III.B, in the presence of multiple modes, physics 
changes. 

A fully numerical multi-mode analysis of the spin-FET 
was proposed by Mireles and Kirczenow (Mireles and 
Kirczenow, 2001). They investigate the effect of the 
strength of the RSOI on the spin-transport properties 
of narrow quantum wires of width W. This is a quasi- 
one-dimensional wire, i.e., assumed to be infinitely long 
in the propagation direction (see figure 11). 

The nearest-neighbour tight-binding Hamiltonian the 
RSOI (5) takes the form 


n/tb 

n so 


(y) = (' 


£,m,s 


' ( is v) s 


C(',m,s C i,j,s‘ (J S 2:)s,s' c t' ,j' 




(47) 


where the summation is intended over next-neighbours 
(.. .) Xi y, tso = f a is th e isotropic next-neighbour transfer 
RSOI hopping term (where a is the lattice constant), and 
cj • s the electron creation operator on the site (£,j) with 
spin state s (s =t, !)• 

The wire is divided in three main regions. In two of 
these [I and III in figure 11], which are near the ferromag¬ 
netic source and drain, the parameter tso is set to zero. 


In the middle region (II) the RSOI is finite (tso 7^ 0). 


The spin-dependent transport problem is solved 
numerically through the use of the spin-dependent 

































16 


Lippman-Schwinger equation 4 


I*) = \®) + G 0 (E)Ug o \y), (48) 


where |4>) is the unperturbed wave function, i.e. an 
eigenstate of the kinetic lattice Hamiltonian Ho with¬ 
out RSOI, and Gq(E) = (E + ie — "Ho) -1 is the Green’s 
functions for the system in the absence of any kind of 
scattering. 

A simple criterion for distinguishing the cases of weak 
and strong RSOI is introduced in Ref. (Mireles and Kir- 
czenow, 2001). The contribution of the mixing of the 
spin sub-bands should be negligible as long the sub-bands 
spacing A% = E ^ — E° is much larger than the sub¬ 
band intermixing energy 


m,s ') 

E° - E° 


< 1 


(49) 


where | (f> n ,s) are the unperturbed electron wavefunctions. 
However, if the confinement energy and/or the RSOI are 
of the same order as the energy shift introduced by the 
sub-band mixing contribution, then the above condition 
is about one or larger. In this case one can introduce 


/?so 


(w) + a ^ F 


PsOi 


(50) 


where /3so = and &f is the Fermi wave number. The 
critical value /3g 0 defines a weak RSOI regime if /3go < 
/3g 0 and a strong coupling regime if f3so > /?so (Mireles 
and Kirczenow, 2001). 

In figure 12 the behaviour of the spin-resolved conduc¬ 
tance as function of the RSOI parameter tgo is shown. 
In figure 12(a) the incident Fermi energy is fixed to 0.5 
(k ~ 0.7a _1 ) and W = 6a = 60 nm, which gives a crit¬ 
ical value /3g 0 = 0.22. This value of /3so separates the 
sinusoidal behaviour of G' 4 for /3go < 0.22 from its be¬ 
haviour for /3go > 0.22 where the confinement energy is 
of the order of the sub-band mixing energy. The effect is 
clearer for a wider wire (W = 120 nm) [see figure 12(c)] 
for which the critical value of /3go is 0.07. To show that 
the non-sinusoidal behaviour is due mainly to the sub¬ 
band mixing, in the figure 12(b) and 12(d) we report 
the spin-conductance as function of tgo with the same 
parameter of figure 12(a) and 12(c) respectively but in 
the unphysical situation of a y ^ 0 and a x = 0, respec¬ 
tively. It is evident that the sinusoidal behaviour is recov¬ 
ered. The modification of spin-FET due to the inclusions 
of many transversal modes has been also investigated in 
Ref. [Jeong and Lee, 2006). 


4 Nowadays a lot of numerical computation, also for systems 
with RSOI, is performed via recursive Green’s function method; 
among the free codes one of the most used is KWANT (Groth 
et al., 2014). 





FIG. 12 Spin-orbit coupling strength dependence of the bal¬ 
listic spin conductance; solid line is G 4 , dashed line G 4 : (a) 
Narrow wire of IT = 6a and uniform RSOI ( a x = a y = 
2atso)- (b) Same as in (a) but with a x = 0 and a y = 2atso; 
perfect oscillations are seen for all tso- (c) Same as in (a) 
with W = 12a. (d) Modulation for W = 12a, with a x = 0 
and a y = 2atso- The sub-band mixing clearly changes the 
otherwise perfectly sinusoidal spin-conductance modulation. 
Courtesy of Ref. [Mireles and Kirczenow, 2001]. 


So far, several obstacles have been found on the way of 
the realisation of the spin-FET (Datta and Das, 1990). 
The main one is related to the injection of spin polarised 
currents. For example, it has been shown that in dif¬ 
fusive transport regime, for typical ferromagnets only a 
current with a small polarisation can be injected into a 
semiconductor 2DEG with long spin-flip length even if 
the conductivity of semiconductor and ferromagnet are 
equal (Schmidt et al., 2000). This situation is dramat¬ 
ically exacerbated when ferromagnetic metals are used; 
in this case the spin polarisation in the semiconductor is 
negligible. 

A possible solution to this problem was proposed by 
employing dilute magnetic semiconductor (Ohno, 1998) 
as source and drain contacts. In these systems a few 
percent of the cations in the III-V or II-VI semiconduc¬ 
tors compounds are randomly substituted by magnetic 
ions, usually Mn, which have local magnetic moments. 
The effective coupling between these local moments is 
mediated by free carriers in the host semiconductor com¬ 
pound (holes for p-doped materials and electrons for n- 
doped one) and can lead to ferromagnetic long-range or¬ 
der. Curie temperatures T c close to 100 K have been 
found in bulk (Ga,Mn)As systems (Ohno, 1998). 
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FIG. 13 (Colours online) (a) Gate voltage modulation of spin 
FETs having different channel lengths, with T = 1.8 K and I 
= 1 mA. The symbols indicate experimental data. The solid 
lines are the fits obtained from V = A cos(2m*aLft _1 + ip). 
Data are offset for clarity. Baseline voltages are 1.032 mV 
and 0.715 mV for L = 1.25 /im and 1.65 /im, respectively; (b) 
Temperature dependence of oscillatory conductance with L = 
1.25 fi m and 1=1 mA. As temperature increases, the mean 
free path decreases and transport characteristics change from 
ballistic to diffusive. Courtesy of Ref. [Koo et al., 2009]. 


Using the properties of the dilute magnetic semicon¬ 
ductor, all-semiconductor spin-FET have been proposed 
in which the conducting channel is provided by a two- 
dimensional hole gas (Pala et al., 2004). 

The first experimental realisation of the spin-FET has 
been done in 2009 by the Johnson’s group (Koo et al., 
2009). Instead of using the configuration initially sug¬ 
gested by Datta and Das (Datta and Das, 1990), in this 
experiment a nonlocal measurement scheme has been 
used (Johnson and Silsbee, 1988; Jedema et al., 2001, 
2002; Lou et al., 2007). This allows to perform a bet¬ 
ter measurement of the charge and spin signals. The 
quantum wire was realised in a InAs heterostructure 
with strong RSOI with the ferromagnetic contacts re¬ 
alised with NigiFeig permalloy on the top. Figure 13(a) 
shows the oscillation of the output voltage of the spin- 
FET as a function of the gate voltage that is changing 
the strength of the RSOI for two different length of the 
distance between the two ferromagnetic contacts. As re¬ 
ported in equation (46) the phase shift between the two 
spin channels is proportional to the distance between the 
two ferromagnets. Therefore, for the longer case (red- 
solid data) we observe a shift of the length on the half 
period of oscillations. Figure 13(b) shows the oscillation 
for a fixed distance between the ferromagnetic contacts 
but at different temperatures. We observe a signature 
of coherent oscillations up to 40 K. At higher tempera¬ 
tures inelastic scattering become more pronounced and 
coherent effects are washed out (Koo et al, 2009). An¬ 
other experimental realization was also lately reported in 
Ref. [Chuang et al, 2014]. 


V. INTERFERENCE EFFECTS AND BERRY PHASE 


A. Rashba interaction as SU(2) gauge field: application to 
quantum networks 


In this section we explore RSOI as a SU(2) non-Abelian 
gauge field. We start by recasting the RSOI in the Hamil¬ 
tonian (8) as a SU(2) vector potential: 


Ho= 2^ 


m a 


2 m* 


P+ X s) 


{p + Asf - 


m a 


(51) 


comparing the two lines of (51) we recognise A s = 
x s). In this form the Rashba Hamiltonian re¬ 
sembles that of a particle in a magnetic field, thus al¬ 
lowing for a straightforward connection with the physics 
of the AB effect (Aharonov and Bohm, 1959). In the 
case of a SU(2) vector potential, it is known as AC ef¬ 
fect (Aharonov and Casher, 1984). In analogy to AB 
effect, we can introduce a phase field that reads: 


0 AC = 


27T 

00 


A s ■ dr = 27r 


0SOI 

00 


(52) 


where 0soi is the flux associated to the RSOI effective 
field and 0o = — the flux quantum. When travelling 
through a closed path an electron gains a non-Abelian 
phase due to the presence of the RSOI. This extra phase 
can give rise to interference phenomena. By contrast 
to the phase gained in a perpendicular magnetic field 

- that depends only on the area enclosed by the par¬ 
ticle path— here the actual path covered by the elec¬ 
tron plays an important role. In the standard AB set-up 
the magnetic field can be tuned in order to move from a 
full destructive to a full constructive interference. In the 
case of AC the role of the magnetic field is played by the 
RSOI that can be modulated modifying a, i.e. by gating 
the heterostructure. The AC effect has been observed in 
semiconductor heterostructure (Bergsten et al, 2006), in 
HgTe rings (Konig et al, 2006) and TI interferometric 
structures (Qu et al, 2011) (c.f. figure 4). The inter¬ 
play between the AB and the AC effect, and the mutual 
effect of an Abelin and non-Abelin gauge field has been 
investigated in Ref. [Nagasawa et al, 2012, 2013]. 

In order to understand the fundamental difference be¬ 
tween the AB and the AC effect, we propose here a very 
simple gedankenexperiment for a square interferometer. 
Let lZ pq be the phase gained by the wave function 'F(r) 

— travelling from a point p to a point q. It reads 


74 u(2 ) =exp{-i/ s ■ (z x dZ)fcso)j 


(53a) 




(53b) 
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FIG. 14 (Colours online) (a): Scheme of the interferometer 
for addressing the difference between an Abelian and a non- 
Abelian phase, (b): Normalised transmission probability T oc 
Trfrr 1 ] for the interferometer in the left panel in the case of 
the action of a SU(2) gauge field (53a). 


depending on whether we are considering a U(l) or SU(2) 
gauge field, where A is the vector potential associated to 
a magnetic field B. We now consider the closed path in 
figure 14(a) and explicitly consider the phases gained by 
the electron travelling along each of its arm. The four 
arms of the interferometer have the same length i pq = L 
(pq £ {AB,AC,BD,CD}), the angle at its base is 9, and 
its area is S. In order to obtain destructive interference 
the sum of the phases gained along consecutive arms has 
to be zero: 

r = <D-<B+^D-<c = 0 (54) 

with /3 e {SU(2),U(1)}. 

By using the Gauss’ theorem and defining the mag¬ 
netic held flux as <f> = |S|5 = (V x A)S, we obtain that 
in the case of a U(l) phase the destructive interference is 
verified when the condition ^ | is met for every the 

value of the angle 6 . Here, the condition for the destruc¬ 
tive interference is straightforward because the phases 
(53b) are c-numbers. 

In the case of the non-Abelian phase held the phases 
(53a) are linear combination of Pauli matrices. There¬ 
fore a solution for the condition (54) is obtained by look¬ 
ing at the transmission probability — proportional (to 
lowest order) to TrpTl]. This quantity is shown in the 
hgure 14(b), we can see that the condition of zero trans¬ 
mission probability is obtained if 9 = \ and k^oL = 

To summarise, while in the case of a U(l) gauge held we 
can obtain destructive interference for each geometrical 
realisation of the interferometer [see hgure 14(a)], in the 
SU(2) case a specific geometric arrangement of the loop 
is required. 

These effects can have interesting implications in the 
transport properties of nontrivial quantum networks. In 
the following, we shall focus on a very special type of lat¬ 
tice structure: the T 3 lattice [c.f. hgure 15(b)], whose unit 
cell contains three inequivalent sites (Sutherland, 1986). 
Its spectrum is composed of two dispersive bands that are 
equivalent to the one of the honeycomb lattice and a hat 
band at zero energy that is due to the presence of lattice 
sites with uneven coordination number (Bercioux et al, 
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FIG. 15 (Colours online) (a) Sketch of diamond chain, this is 
the one-dimensional version of the T 3 lattice, (b): Sketch of 
the 7i lattice, it is a Bravais lattice with a unit cell containing 
three inequivalent sites: A, B and H, the latter has coordi¬ 
nation number higher that the other two (Sutherland, 1986). 

(c) : averaged conductance for the diamond chain as a func¬ 
tion of the magnetic flux ratio (lower oxaxis, blue solid-line) 
and as a function of the dimensionless RSOI (upper oxaxis, 
red-dashed line) (Courtesy of Ref. [Bercioux et al, 2005b]). 

(d) : averaged conductance normalised to the number of leads 
for the 75 lattice as a function of the magnetic flux ratio 
(lower er-axis, blue solid-line) and as a function of the dimen¬ 
sionless RSOI (upper axaxis, red-dashed line) (Courtesy of 
Ref. [Bercioux et al, 2005b] ). In (c) and (d) B and E are the 
applied magnetic field and the electric field, respectively. 


2009, 2011). In the presence of a perpendicular magnetic 
field — when the ratio ^ = x — the spectrum reduces 
to three flat bands (Vidal et al., 1998). This localisa¬ 
tion effect is a consequence of the AB effect taking place 
within each of the plaquette contained in the T 3 lattice 
[in yellow in figure 15(b)]. The same effect is observed in 
its quasi-ID analog, the diamond chain [c.f. figure 15(a)] 
for the Jjj- = j (Vidal et al., 2000a). 

Can the same physics be driven also by the RSOI? The 
question has been addressed by one of the authors using 
the method of quantum networks (Naud et al., 2001). 
Full localisation due RSOI is possible only in the dia¬ 
mond chain (Bercioux et al., 2004) [figure 15(c)] and is 
forbidden, for the arguments illustrated above, in the case 
of the T 3 lattice (Bercioux et al., 2005b) — figure 15(d). 
In figure 15(c) and 15(d) we show the conductance of 
the diamond lattice and the T 3 lattice, respectively, as 
a function of the magnetic flux and RSOI. We observe 
that in the case of the diamond chain 15(c), both effects 
induce a complete localisation equivalent to zero conduc¬ 
tance. However, in the case of the T 3 lattice 15(d), only 
the magnetic field can induce (almost) complete locali¬ 
sation. The residual conductance comes from the con- 
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ducting state along the lattice boundary due to the Hall 
effect. However it can be reduced to zero by injecting 
electrons through the centre of the lattice (Vidal et al. , 
2000b). Signatures of this AB caging effect have been 
observed also experimentally (Naud et al., 2001; Abilio 
et al., 1999). The diamond chain lattice model with RSOI 
and AB field has been also proposed as a system for pro¬ 
ducing a spin filter (Hatano et al., 2007; Aharony et al., 
2008, 2011; Matityahu et al., 2013a, b). 


B. Interference in quantum rings 

Mesoscopic quantum rings allow to have direct access 
to the phase of the electron wavefunction, when their 
size is smaller than the coherence length. Interference 
effects have been observed in metal quantum rings many 
years ago (Washburn and Webb, 1992). Since electrons 
are spinful particles, the spin part of the wavefunction is 
influenced by the magnetic field via the Zeeman term in 
the Hamiltonian. A more subtle effect arises when there 
is a magnetic field non-orthogonal to the plane of the or¬ 
biting particle ( e.g. the effective magnetic field due to 
the RSOI) because, as a consequence of the orbital mo¬ 
tion, its spin dynamics is instantaneously governed by a 
time-dependent Hamiltonian (Anandan, 2002; Saarikoski 
et al., 2014). This time dependence ends up in an ex¬ 
tra phase acquired by the particle wavefunction which 
is named after Berry [Berry, 1984 and Aharonov and 
Anandan, 1987], who put in foreground its topological 
properties when the orbits are closed. 

During the last years, the effects of RSOI on the AB 
oscillations have been observed in semiconductor based 
quantum rings by several groups (Yau et al., 2002; Nitta 
et al., 2003; Meijer et al, 2004; Morpurgo et al, 1998). 
As said before, in the presence of both orthogonal mag¬ 
netic field and RSOI, the total effective momentum de¬ 
pendent magnetic field is tilted with respect to the ver¬ 
tical direction. The resulting Berry phase influences the 
interference pattern. 

In Refs. [Capozza et al., 2005 and Lucignano et al., 
2007] one of the authors studied the conductance and the 
spin transport in a quantum ring in the presence of RSOI 
and magnetic field, accounting also for dephasing at the 
contacts. The AB resonance in the Fourier transform of 
the magneto-conductance displays satellite peaks due to 
the RSOI (see figure 16) that have been experimentally 
observed (Habib et al., 2007). During the last years, sev¬ 
eral theoretical techniques have been employed to study 
quantum rings. In Refs. [Loss et al., 1990 and Engel and 
Loss, 2000], an imaginary time path integral approach is 
developed to study the conductance of a strictly ID quan¬ 
tum ring, and its conductance fluctuations in the diffusive 
limit. In Ref. [Tserkovnyak and Brataas, 2007] a real time 
path integral approach is applied in the limit of negligible 
Zeeman splitting. Several research articles discussed the 


conductance properties and the spin dependent transport 
of quantum rings in the ID ballistic limit, by means of 
a spin dependent scattering matrix approach (Frustaglia 
and Richter, 2004; Bercioux et al., 2005a; Molnar et al., 
2004; Citro et al, 2006; Ramaglia et al, 2006). In the ab¬ 
sence of the magnetic flux, the conductance shows quasi- 
periodic oscillations in the RSOI strength, which can be 
modified by switching the magnetic field on. Numerical 
calculations (Frustaglia and Richter, 2004; Sourna and 
Nikolic, 2004; Wu and Cao, 2006) have shown that in the 
2D case there are only quantitative modifications of the 
ID results that do not qualitatively affect the physics. 
Thus, in the following, we will just focus on the ID limit. 

The model Hamiltonian, describing a half integer 
spin particle, in an orthogonal magnetic field and with 
RSOI (Meijer et al., 2002), reads 


n\p,r,S\ = 


1 

2 m 


(p + — uj c S z + Hrsoi (55) 

(p+ 6 ^ 

V c 


c 

2 a r 

'HrSOI = ^2 [ ZX 


where Ao(r) = ^-(—y,x,0) is the vector potential gener¬ 
ating the uniform field B, normal to the ring surface, and 
w c = 9 2 eB is the cyclotron frequency. We fix the vector 
potential in the symmetric gauge. We present, here, fully 
general results as they depend on the ratio R which 
can be tuned by acting on a. We assume a single channel 
ring as a ID circle of radius R, connected to two leads. 
Accordingly, the position of the particle on the ring is 
parametrised by the angle <p. The vector potential has 
just the azimuthal component A v = 5 )^, where <j> is the 
magnetic flux threading the ring. 



FIG. 16 (Colours online) a): Sketch of a quantum ring. The 
magnetic field B drives the AB effect, while the electric field 
can tune the RSOI. b) Fourier transform of the magneto con¬ 
ductance. By increasing ROI two satellite peaks rise close to 
the AB peak. They can be attributed to the Berry phase 
acquired by the electron during its motion. Courtesy of 
Ref. [Capozza et al, 2005]. 
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In order to study the conduction properties of the ring, 
we need the propagation amplitude for an electron at 
energy Eq entering the ring with spin polarisation sq and 
leaving it with spin polarisation sf . This is given by 

A(sf,s 0 \E 0 ) = / (r f ,Sf,tf\r 0 ,s 0 ,t 0 ), (56) 

J o T o 

where (ry, s/, f/|ro, sg, to) is the amplitude for a parti¬ 
cle entering the ring at the point r o and at the time t o 
with spin polarisation sq to exit at the point rf at the 
time tf with spin polarisation Sf. Here, tq = r ”)) is the 
time scale for the quantum motion. In order to compute 
(rf,Sf,tf\ro,so,to), we adopt a path integral represen¬ 
tation for the orbital part of the amplitude. Since we 
parametrise the orbital motion of the particle in terms 
of the angle ip, we provide the appropriate Lagrangian, 
£ or b, as a function of (p, dp. It reads 

A>rb [<p(t),<p(t),tr\ = ™R 2 ip 2 (t)- j~Jtp(t) 

a 2 m h 2 

+ ~2fR + &mR? ■ (5?) 

The last two contributions in equation (57) are constants 
that come from the RSOI term, and the Arthurs term, 
which is required when a path integration is performed 
in cylindrical coordinates. Since both contributions are 
constant, they can be lumped into the incoming energy 
Eq ~ Ep and therefore they will be omitted henceforth. 

By taking into account the spin degree of freedom, as 
well, we represent the propagation amplitude as 

( rf,Sf,tf\r 0 ,s 0 ,t 0 ) = (r 7 ,s/|e -ljt o dl H \r 0 ,s 0 ,t 0 ) 

n(f(t f)=ip f t 

= Vpe-'to dt [^ov 2 -H ( Sf \U spin (t f ,to)\s 0 }, 

Jtp(t 0 )=tpo 

(58) 


where q = 


^Apin (If : +i) — T exp 


rtf 


df Jospin) 7 ') 
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(59) 


is the full spin propagator, T is the time ordering oper¬ 
ator, and H spin (t) is the spin Hamiltonian that reads 


H 


spin 


(t) 


hf uj c 

2 ^ 2y(pe i vW -ui c 


(60) 


with 7 = . 

For a ring device, at each lead one has to take into 
account three possible scattering processes, consistently 
with the conservation of the total current. This is de¬ 
scribed in terms of a unitary S —matrix that, when the 
two arms are symmetric, is given by 


5 = 


( — 1(1 + r) 1(1 -r) \j\(l~r 2 ) ^ 

i(!-r) -1(1 + 0 \J\(l-r 2 ) 

\ \]W-r 2 ) r 


(61) 


The numerical labelling of the 5-matrix elements refer¬ 
ring to the three terminals of each contact fork, are ex¬ 
plained in figure 17(la). Assuming, for simplicity, that 
the scattering matrix is the same for both leads, equa¬ 
tion (61) will hold both at the left-hand lead, and at the 
right-hand lead of the ring. 



FIG. 17 (Colours online) First and second order paths in¬ 
cluded in the calculation of the transmission amplitude across 
the ring, from left to right, including forward scattering only. 
Numbers 1, 2, 3(T, 2', 3') in figure 17a refer to the labelling of 
the terminals in Eq.(61). Courtesy of Ref. [Lucignano et al . , 
2007], 

In particular, £ 3,3 = r is the reflection amplitude for 
a wave coming from the left lead, 5 l( 2 ),i( 2 ) = — ^(1 + r) 
the reflection amplitude for a wave incoming from the up¬ 
per/lower arm, 5 l( 2 ), 2 (i) = |(1 — r) is the transmission 
amplitude from the upper (lower) to the lower (upper) 

arm and ^ 3 ( 2 )^ = 5 3 ,i( 2 ) = \j\( 1 — f 2 ) i s the transmis¬ 
sion amplitude from the upper/lower arm to outside of 
the ring. As the ring is assumed to be symmetric, the 
same scattering matrix applies to the right lead where 
we indicate with primed numbers the three terminals as 
shown in figure 17. In figure 17 we show the simplest 
paths of the electrons in the ring including only forward 
scattering at the contacts. More involuted paths arise if 
we account also for backscattering processes in which the 
electron can get backscattered within the same ring’s arm 
from which it is coming. For instance, the paths (2/) and 
( 2h ), as well as ( 2 <?) and ( 2 i) in figure 18, include looping 
in opposite directions around the ring. Interference be¬ 
tween clockwise and counterclockwise windings leads to 
WL corrections. We denote these corresponding paths 
— including also (2c) and (2d) — as “reversed paths”. In 
our approach, all order paths are numerically generated 
up to the convergency and the S— matrix (61) is imple¬ 
mented in the numerical algorithm. Within Landauer’s 
approach, the conductance Q is given by 

£ = \As;s'\e f )\ 2 . (62) 

s,s' 

Here we will consider the dependence on the external 
magnetic field ^ and on the RSOI strength ksoR both 
in the absence and in the presence of dephasing at the 
contacts. To make the model more realistic, we allow for 
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FIG. 18 Second order paths of the transmission amplitude 
from left to right including backscattering at the leads. Paths 
(2/) and (2 h), as well as (2 g ) and (2 i) contribute to the weak 
localization corrections. Courtesy of Ref. [Lucignano et al, 
2007], 

higher order looping of the electron within the ring. In 
Ref. [Capozza et al. , 2005], only the paths of the kind 
of figures 17(2a) and 17(2b) were included. Following 
Ref. [Lucignano et al. , 2007] we consider here also the 
paths of the kind of figure 18 in which the electron can 
be backscattered into the ring. We use here r = 0 in 
the scattering matrix between the arms and the leads, 
which means that no back-reflection in the incoming lead 
is present. 

The dephasing due to diffusiveness in the contacts is 
accounted by adding a random phase z £ (—£,£) for 
each scattering at the leads. In figure 19, we report the 
conductance Q as a function of ^ , with ksoR = 0 (left 
column) or as a function of ksoR with = 0 (right 
column). These are averaged over N = 1000 realisations 
of dephasing, and plotted for increasing window of phase 
randomness (£ = 7r/3,7T,27r from top to bottom). The 
black curves refer to ideal contacts (including only the 
paths of figure 17) while the red curves refer to realistic 
contacts (including also the paths of figure 18). 

The ring is rather insensitive to small dephasing at the 
contacts however, by increasing the amount of dephasing 
(middle and bottom left panels in figure 19) we find that 
the sensitiveness is larger in the case of realistic contacts. 
This is due to the fact that for realistic coupling, the 
electrons in the ring can experience higher order paths, 
since it scatters with the leads many times. 

In the right panel of figure 19, we plot the dc- 
conductance as a function of ksoR at (f> = 0 for both ideal 
contacts and realistic contacts (black and red lines in 
each box), with an increasing phase randomisation (boxes 
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FIG. 19 (Colours online)Conductance as a function of <j>/<j>o 
(left panels) and ksoR (right panels) for ideal (black curves) 
and realistic (red curves) contacts. An increasing amount 
of dephasing at the contacts is also included: from top to 
bottom: £ = n/3,n,2n. Courtesy of Ref. [Lucignano et al, 
2007], 

from top to bottom with £ = ^,7r, 27 t), averaged over 
N = 1000 disorder realisations. In the case of ideal con¬ 
tacts and little dephasing (top right panel black curve), 
we observe again quasi periodic oscillation of the con¬ 
ductance reproducing, the localisation conditions at the 
expected values of ksoR (Frustaglia and Richter, 2004; 
Bercioux et al., 2005a; Capozza et al., 2005; Lucignano 
et al., 2007; Molnar et al, 2004). When including higher 
order processes, interference effects give rise to a slightly 
different pattern. In the case of realistic contacts, we note 
that the conductance of the ring is seriously affected by 
dephasing. Indeed, large dephasing gives rise to random 
oscillations that are not averaged out — thus washing out 
the conductance oscillations. The effect takes place for 
£ ~ 7 r when time reversed paths are included, in contrast 
to £ ~ 27 t when the time reversed paths are absent. As 
regular magneto-conductance oscillations are experimen¬ 
tally observed (Yau et al., 2002; Nitta et al., 2003; Meijer 
et al, 2004; Morpurgo et al., 1998) with little percentage 
of contrast between maxima and minima, we conclude 
that, in real samples, dephasing is ubiquitous. 

VI. RELATED PROBLEMS 
A. Spin-Hall effect 

The SHE is a phenomenon, associated to SOI — it can 
be used to electrically generate or detect spin currents 
in non-magnetic systems. This effect has been observed 
both in metallic and semiconductor systems (for a review 
see Ref. [Jungwirth et al., 2012 and Vignale, 2010]). In 
this review, we will mainly focus on semiconductors, in 
which RSOI plays a major role. The SHE predicts that 
an unpolarised electric current can generate a transverse 
spin current, whose spin is perpendicular to the plane 
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of the two currents. It was predicted in two pioneering 
articles by Dyakonov and Perel in 1971 [D’yakonov and 
Perel, 1971a, b]. It is a consequence of the Mott scatter¬ 
ing of electrons on unpolarised impurities, which results 
in spatial separation of electrons with opposite spins. It is 
closely related to the anomalous Hall effect — Hall effect 
in ferromagnetic materials — originally observed by Hall 
himself [Hall, 1881] and later explained by Karplus and 
Luttinger [Karplus and Luttinger, 1954] and Nozieres and 
Lewiner [Nozieres and Lewiner, 1973]. It does not require 
magnetic field nor magnetism, in other words it does not 
require broken time reversal symmetry (for a review see 
[Nagaosa et al. , 2010]). The SHE was hardly investigated 
till Hirsch [Hirsch, 1999] and Zhang [Zhang, 2000] pro¬ 
posed such phenomenon to the attention of the spintron- 
ics community. Three main mechanisms have been pro¬ 
posed to contribute to the SHE: namely the spin depen¬ 
dent band structure of the material (the so called intrin¬ 
sic mechanism) and the impurity scattering mechanisms, 
i.e. the “skew scattering" and the “side jump", conven¬ 
tionally addressed as extrinsic mechanisms. The extrinsic 
mechanisms are the same responsible of the anomalous 
Hall effect (Nagaosa et al., 2010): the skew scattering 
originates from the different scattering angle due to a 
spin-orbit impurity, depending on the electron spin and 
angular momentum; while the side jump reflects the shift 
in the trajectory of the backscattered electrons originat¬ 
ing from the anomalous velocity operator in SOI systems 
(Crepieux and Bruno, 2001). The intrinsic mechanism 
was proposed by Murakami et al. [Murakami et al., 2003] 
and by Sinova et al. [Sinova et al., 2004], it resorts on 
the DSOI and RSOI in 2D III-V GaAs quantum wells. 
This mechanism does not explicitly depend on impurities, 
however these cannot be neglected as they are absolutely 
essential to the establishment of the steady state current. 
In principle, the spin Hall current can be calculated using 
the Kubo formula. However, inclusion of vertex correc¬ 
tions due to impurities, in the linear SOI model leads to a 
vanishing spin Hall conductivity (Raimondi and Schwab, 
2005; Dimitrova, 2005) in striking contrast to earlier cal¬ 
culations not including vertex corrections (Sinova et al, 
2004). More complicated models (explicitly carried out 
for transition metals) (Tanaka et al., 2008) overtake this 
problem, giving rise to a finite spin Hall conductivity in 
reasonable agreement with the experiments. Some pro¬ 
posals require the presence of a magnetic field to give 
rise to a finite spin Hall current even in the presence 
of linear SOI and disorder (Lucignano et al., 2008; En¬ 
gel et al, 2007; Lin et al, 2006; Milletarf et al., 2008; 
Raimondi et al., 2006; Gorini et al., 2008). Among the 
first successful experiments we here highlight the work 
by Wunderlich et al. [Wunderlich et al., 2005c] that used 
coplanar p-n diodes to detect circularly polarised electro¬ 
luminescence at opposite edges of the spin Hall bar, and 
the work by Kato et al. [Kato et al., 2004] that employed 
a magneto-optical Kerr microscope to scan the spin po¬ 


larisation across the Hall bar. The former ascribed their 
finding to the intrinsic mechanism whereas the latter to 
the extrinsic one. Since that time several interesting ex¬ 
periments performing optical measurements for the spin 
detection in the intrinsic and extrinsic SHE have been 
performed (Matsuzaka et al., 2009; Nomura et al., 2005; 
Sih et al., 2005b, 2006b; Stern et al, 2006, 2007, 2008; 
Zhao et al, 2006; Okamoto et al., 2014). A detailed de¬ 
scription of this physics is out of the purpose of this re¬ 
views and we refer the readers to more specific review ar¬ 
ticles on the SHE for more information (Jungwirth et al, 
2012 ; Vignale, 2010). 


B. Diffusive limit: weak anti-localization 


In the previous sections we have mainly described the 
physics of ballistic mesoscopic systems, i. e. solid state de¬ 
vices in which the lateral size of the sample L is smaller 
than the coherence length L^ and of the elastic mean free 
path t e . In this regime, impurity scattering can be safely 
disregarded, and the ballistic motion of electrons well 
describes the physical scenario. However, in the meso¬ 
scopic limit, it may happen that the lateral size of the 
sample L, while being still smaller (or of the order of) 
than the coherent length of the material L < L</,, it is 
larger than the elastic mean free path: L > l e . Impurity 
scattering cannot be neglected and the electrons scatter 
in the sample undergoing a quantum diffusive dynamics. 
In a 2DEG at low temperature, quantum corrections to 
the conductance are mostly due to interference effects 
between two electrons scattered by the same impurity. 
These propagate along the same closed trajectory but in 
opposite directions — usually known as Cooperon correc¬ 
tion. This process has its more spectacular manifestation 
in the so called WL (see Ref. [Akkermans and Montam- 
baux, 2007; Altshuler et al., 1980; and Kramer, B. and 
MacKinnon, A., 1993]): it increases the effective scat¬ 
tering cross-section and therefore leads to a suppression 
of the conductivity. In the presence of a magnetic field 
B perpendicular to the electron plane, the two waves 
acquire a phase difference A ip = 2^-, where 4>b is the 
magnetic flux through the area enclosed by the electron 
trajectory, therefore quantum interference effects are sup¬ 
pressed and the conductance increases to the classical 
limit. Such a modification in conductivity, as a function 
of the magnetic field, is known as positive magneto con¬ 
ductivity (sometimes called negative magnetoresistance). 
This quantum correction to the resistivity can be explic¬ 
itly calculated evaluating the Cooperon diagram. For a 
2DEG in perpendicular magnetic field 


Aa(B) oc 


/*oo 

/ df Z c (t,B) 

J o 



(63) 


where Z c (t) = ^ sinh 1 ( 47r ^,* B ), where and r e are 
phase coherence time and electron relaxation timeand D 





23 


is the 2D diffusion coefficient, respectively. In the limit 
of very weak magnetic field (B <C s ^ 2 ) the magnetore- 

sistance results to be: 

Ai?(B) oc, (64) 

a Q 

thus negative, and the resistivity correction can be shown 
to be 


Aa(B) — Acr(O) oc 



(65) 


where B$ = is a characteristic magnetic field corre- 

sponding to a flux quantum through an area of the order 
of nL Measuring the negative magneto-resistance by 
applying a weak magnetic field is a very elegant way to 
probe weak localization correction to the conductance. 

In the presence of SOI the WL correction changes its 
sign and the magnetoresistance results to be positive. In¬ 
deed, Cooperon correction is multiplied by a factor 


(Qsoiit)) 



t 

T soi 



where tsoi sets the time scale of the SOI. For strong 
SOI, the exponential factor can be neglected at any time 
t, thus the factor (Qsoi) changes its sign and tends to 
— It results in a change of the sign of conventional 
weak localization correction. This phenomenon, theo¬ 
rised in Ref. (Hikami et ai, 1980), is called WAL. It was 
first observed in metal films in the presence of spin-orbit 
impurities (Bergman, 1984, 1982) and has been inten¬ 
sively studied in the recent past, both in 2DEGs (Miller 
et al., 2003; Knap et al, 1996; Koga et ai, 2002) and in 
large quantum dots (Zumbiihl et al., 2002; Aleiner and 
Fal’ko, 2001) made out of III-V semiconductors in the 
presence of RSOI [c.f. for instance figure 1(b)]. Recently, 
weak WAL measurements have been used to quantify the 
RSOI strength in oxide interfaces (Caviglia et ai, 2010) 
and TIs (Bardarson and Moore, 2013). 


C. Quantum wires with RSOI, superconductivity and 
magnetic field: The quest for Majorana Fermions 

In his seminal article on the symmetric theory of elec¬ 
tron and positron (Majorana, 1937), Ettore Majorana 
predicted the existence of selfadjoint fermions, i.e. half 
integer spin particles being their own antiparticle, as real 
solutions of Dirac equation. Since the first Majorana’s 
prediction, the quest for Majorana Fermions (MFs) as 
elementary particles, in the high energy physics, is still 
open. Experiments on the neutrinoless double beta de¬ 
cay should unveil whether neutrinos are MFs or not, but 
to the date there is no answer to this question. 


On the contrary, electrons in solid, under particular 
circumstances could recombine with holes — their rela¬ 
tive anti-particles in the solid state language — in order 
to form self-adjoint excitations, i.e. Majorana quasipar¬ 
ticles. 

Thus MFs (Kitaev, 2001, 2006) can appear as quasi¬ 
particle excitations in solid state systems. For instance, 
they are expected to show up as boundary states of the 
Kitaev model: a toy-model Hamiltonian describing ID 
spineless p -wave superconductors. 

Despite the fact that “solid state” MFs are not elemen¬ 
tary particles, there is strong excitement in the condensed 
matter community, and many efforts are being devoted 
to their theoretical and experimental understanding. Of 
course, finding a novel elementary excitation in solids 
would have a fundamental significance on its own, how¬ 
ever, in the case of MF there is strong interest also for 
its possible applications. Indeed MFs, due to their non- 
Abelian statistics (Read and Green, 2000; Ivanov, 2001; 
Stern et ai, 2004; Alicea et ai, 2011) can be considered 
as building blocks for topological quantum computation 
and for other alternative decoherence free quantum com¬ 
putational schemes (Nayak et al., 2008). 

Several mechanisms and devices have been proposed 
to isolate and detect MFs in meso/nano scaled devices 
(c.f. Refs. [Alicea, 2012 and Beenakker, 2013] for ex¬ 
haustive reviews on the topic). However, to the date, 
while many efforts are focused on hybrid TI supercon¬ 
ductor heterostructures (Veldhorst et al., 2012; Galletti 
et ai, 2014a, b), (one of) the most promising device re¬ 
sorts on the use of III-V semiconducting nanowires with 
strong RSOI, in the proximity with conventional super¬ 
conductors and in the presence of a magnetic field along 
the wire. This mechanism, early proposed in Refs. (Oreg 
et al., 2010; Lutchyn et ai, 2010) has been later experi¬ 
mentally explored in Refs. [Das et ai, 2012; Mourik et ai, 
2012; and Churchill et al, 2013]. While strong hints of 
MF physics have been highlighted in Ref. [Mourik et ai, 
2012], other recent articles [Churchill et ai, 2013; Lin¬ 
der et al, 2010; Lee et al . , 2012; Liu et ai, 2012; Kells 
et ai, 2012; and Sun and Shah, 2015] offer alternative in¬ 
terpretation of the observed zero bias anomaly in terms 
of Kondo physics. Therefore an unambiguous “smoking- 
gun” experimental proof of the presence of MFs still has 
to come. 

In the following, we will describe the theoretical back¬ 
ground supporting the physics of MFs in quasi-ID nano 
wires with RSOI, magnetic field and proximity induced 
superconductivity. After that, we briefly review a recent 
experimental result. 

For sake of simplicity we assume a single active chan¬ 
nel in the nanowire and set h = 1. The nanowire is along 
the x direction, therefore the transverse dynamics, along 
y direction, can be neglected. The Hamiltonian (15) de¬ 
scribing the wire can be simplified as follows: 







"Hmf = E<t<7' fZ ^t( x ) (Z* - as y^ + 9HbBs x - /z) Vv(z) + + h.c.j da; 
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( 66 ) 


where we have also added a Zeeman term due to the 
magnetic field B along the wire axis and the supercon¬ 
ducting pairing potential A. The operator i/>J.(x) creates 
an electron of spin a at the position x. One can notice 
that the Hamiltonian ( 66 ) is not that of a ID spineless 
p -wave superconductor, indeed we can easily recognise 
a local s-wave superconducting pairing and spinful elec¬ 
trons. However, by analysing the action of the various 
interaction terms in ( 66 ), we can show how to obtain 
an optimal choice of the parameter regime that would 
map ( 66 ) into the effective Hamiltonian of a ID spinless 
p -wave superconductor. 

We choose 'F(x) = [u-|-(x),U|(x),^(x), —r>-|-(x)] as a 
Nambu spinor — where u/v are quasi electron/hole wave 
functions of momentum k (De Gennes, 1999). Thus the 
Bogolubov-de Gennes (BdG) Hamiltonian ( 66 ) in the 
'F(x) basis, can be recast in the simple matrix form: 


%BdG ( k ) = 


( k 2 


\2m* 


+ i aksy - /i ) t z + E z s x + A t x (67) 


where the n Pauli matrices act in the Nambu space, 
while Si act in the spin space and we have introduced 
the Zeeman energy E z . In the absence of magnetic field 
and superconductivity, the conventional spin-degenerate 
parabolic band dispersion splits because of RSOI, into 
two parabolas crossing at k = 0 [c.f. figure 5(e)]. The 
Zeeman coupling, splits the two bands at k = 0 giving 
rise to the two bands E + {k) and E_{k) in figure 20(a). 
If the Fermi energy lies between the two bands, and the 
temperature is low enough that the higher band cannot 
be thermally populated, the system is effectively a ID 
spinless system. The spin degree-of-freedom has been 
effectively quenched by the simultaneous action of the 
RSOI and Zeeman splitting. Only a single pseudo spin 
degree-of-freedom “ —" is involved here. In addition we 
have the s-wave superconducting pairing. Projecting the 
full Hamiltonian onto the E_{k) band we can write down 
an effective Hamiltonian 

Up = Y^ (b~ i k )c\c k + A_ (fc)44 + h.c.) ( 68 ) 

k 


where A_(fc) = \ak x A{a 2 k 2 + E ?) -1 ^ 2 is an effective 
superconducting pairing with the desired p -wave symme¬ 
try. Thus, the present problem is isomorphic to the Ma- 
jorana wire considered by Kitaev [Kitaev, 2001]. That 
guarantees the presence of MFs as boundary excitations. 
However in the spirit of a review article, here we will not 
provide an explicit proof of the presence of MFs as bound¬ 
ary excitations of the wire, as it can be inferred by the 
isomorphism to Kitaev model and can also be found in 
Refs. [Oreg et al., 2010 and Lutchyn et al., 2010] but we 


will provide the reader with simple symmetry arguments 
to qualitatively justify this result. 



FIG. 20 (Colours online) (a) Free electrons band splitter by 
the simultaneous application of RSOI and Zeeman splitting, 
(b) Phase diagram of the model Hamiltonian (67). 

Let us start by commenting on the spectrum. By 
squaring twice Hamiltonian (67) we obtain the disper¬ 
sion relation for the two bands 

E 2 ± (k, A) =E 2 + A 2 +e k + (ak) 2 

±2^/{E z AY + (E z £ k ) 2 + (a/c^) 2 

where = k 2 /2m — p. The pairing A plays two crucial 
roles. It opens a gap at the outer wings of the dispersion, 
where the Zeeman field is unimportant, and modifies the 
gap forming near k = 0. The former role eliminates the 
possibility of high-momentum gapless excitations, thus 
leaving only the chiral states near k = 0 as low energy 
excitation. The latter role allows us to tune the topolog¬ 
ical phase transition essential for isolating MFs. One can 
easily verify that the Hamiltonian (67) satisfies particle- 
hole symmetry: 

= Ub<ig(— k) ■ (69) 

Hence, for each eigenstate of positive energy Ubag^ = 
Ei/j it exists a corresponding eigenstate of opposite 
energy "HedG^^) = —E(Aip). A linear combination 7 
. + (Ety>) is of course a selfadjoint Fermionic operator, 
i.e. a MF, however it is not an Hamiltonian eigenvector, 
unlike E = 0 (Chamon et al. , 2010). 

It can be shown that the presence of such zero energy 
selfadjoint eigenstate is indeed related to the sign of the 
gap at k = 0 

E 0 = \E + (k, A) - E-(k, A )| fc=0 = | E z - 4A 2 + p 2 \ 
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FIG. 21 (Colours online) (a): Experimental device. Electrons 
from a normal metal are injected in the InSb nanowire which 
is in contact with NbTiN superconductor, (b): Differential 
conductance traces as a function of the magnetic field. As 
soon as the topological criterion is satisfied, a zero bias peak 
appears. Courtesy of Ref. [Mourik et al, 2012]. 


between the two excitation branches E± (k, A). When the 
Zeeman energy E z closes such gap, we have a quantum 
phase transition between two different physical regimes. 
At E z > a/A 2 " + /i 2 the gap is a Zeeman gap while in 
the opposite case the gap is due to superconductivity. 
Explicit construction of the zero energy excitation in the 
two cases shows that when the gap is magnetic, selfad- 
joint zero energy excitations, i.e. MFs can appear at the 
two boundaries of the wire. The phase diagram of such 
system is shown in figure 20(b). Out of the topological 
regime no MFs are expected, but the interesting physics 
of anomalous Josephson effect sets in (Yokoyama et al ., 
2014; Campagnano et al., 2014). 

In Ref. [Mourik et al., 2012] electrons are injected from 
a normal metal in an InSb nanowire in proximity with 
a NbTiN superconductor [see figure 21(a)] and the dif¬ 
ferential conductance is measured as a function of the 
magnetic field. As proposed in Ref. [Flensberg, 2010] 
electrons from the normal metal scatter against the zero 
energy MF thus revealing a zero bias anomaly in the dif¬ 
ferential conductance. In figure 21(b) we can see that as 
soon as the magnetic field is large enough to satisfy the 
topological criterion, a zero bias peak appears that is in¬ 
terpreted as the resonant scattering of electrons through 
the zero energy Majorana state. Experiments on these 
kind of structures are very challenging: on the one hand 
in order to satisfy the topological criterion, a sizable mag¬ 
netic field is required, on the other hand, the same mag¬ 
netic field can disrupt the very fragile superconductivity 
induced in the semiconducting nanowire by the proxim¬ 
ity effect (Potter and Lee, 2011). Everything has to be 
finely tuned and the temperature has to be quite low, 
usually below lOOmK. That is why many other proposals 
involving high critical temperature superconductors (Lu- 
cignano et al., 2012, 2013) or completely different plat¬ 
forms (Nadj-Perge et al., 2014; Hart et al, 2014) have 
been investigated, too. 


VII. CONCLUSIONS AND OUTLOOK 

In this review, we have inspected some of the most in¬ 
teresting theoretical and experimental results on spin de¬ 
pendent quantum transport in mesoscopic systems in the 
presence of RSOI. During the last decades a significant 
number of theories and experiments have appeared. The 
main goal of this review is to give an historical overview 
of the field so to address the reader directly to a more spe- 
cialistic literature. Rashba spin-orbit interaction turned 
out to be an invaluable tool to have access and to manip¬ 
ulate the electron spin degree of freedom, without break¬ 
ing time reversal symmetry. From the applicative point 
of view, RSOI has been intensively studied to generate 
spin currents in semiconducting electron systems. Sev¬ 
eral mechanisms and devices for producing pure spin cur¬ 
rents have been put forward. Among them, here we have 
focused our attention on quantum spin ratchet and spin 
pumping. 

The physics associated to the RSOI can have very ap¬ 
plicative oriented aspects but can also be of stimulus 
for addressing more fundamental issues of quantum me¬ 
chanics, e.g., quantum interference effects. These have 
been shown to be strongly affected by the RSOI, with 
the novelty of the phase factors of non-Abelian nature. 
We have seen how this additional phase can give rise to 
unexpected phenomena, such as the localisation in net¬ 
works with nontrivial connectivity, or the appearance of 
anomalous peaks in the Fourier transform of the mag¬ 
neto conductance in quantum rings. Recently, new and 
exciting perspectives open, because RSOI is of funda¬ 
mental relevance for the physics of topological insulators 
and of Majorana fermions in hybrid superconducting¬ 
semiconducting heterostructures that are currently at¬ 
tracting a lot of interest not only in the condensed matter 
community. 

Besides charge and spin transport, SOI may have a rel¬ 
evant role also in heat transport which is the subject of 
spin caloritronics (Bauer et al, 2012). Here, the com¬ 
bination of magnetic and non magnetic hybrid structure 
can lead to anomalous behaviour in the Seebeck and/or 
the Peltier coefficients (Gravier et al, 2006; Flipse et al, 
2012). Deviation from standard behavior has been ob¬ 
served also in experiments on SLG (Zuev et al, 2009; 
Wei et al, 2009). Furthermore, it is well known that 
the most efficient media for creating thermoelectric de¬ 
vices are based on materials with a very strong SOI as 
HgTe, PbTe, Bi,Te, and Bi/Sb alloys, Ref. (Delves, 1965; 
Biswas et al, 2012). However, a complete understand¬ 
ing of the role of RSOI on these properties is still under 
investigation (Alomar and Sanchez, 2014; Rameshti and 
Moghaddam, 2015). 

We are sure that in the future this field of research 
will be still thriving. Most probably new electronic and 
spintronic devices will come to the market with function¬ 
ality that are associated to the physics of the RSOI. On 
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the other side, as it happened for the case of the Majo- 
rana quasiparticle research, RSOI could pave the way to 
fundamental research in order to investigate exotic phe¬ 
nomena that could get accessed by quantum simulators 
realised in solid state systems (Feynman, 1982). 
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